WORLD INTELLECTUAL PROPERTY ORGANIZATION 
International Bureau 




PCT 

INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(51) International Patent Classification 6 : 

G01N 33/53, 33/573, 33/574, G06F 
159AM) 



Al 



(11) International Publication Number: WO 96A2187 

(43) International Publication Date: 25 April 1996 (25.04.%) 



(21) International Application Number: PCT/US95/01 379 

(22) International Filing Date: 2 February 1995 (02.02.95) 



(30) Priority Data: 
323,446 



13 October 1994(13.10.94) 



US 



(71) Applicant: HORUS THERAPEUTICS . INC. [US/US]; 2320 
Brighton-Henrietta Town Line Road. Rochester, NY 14623 
(US). 

(72) Inventors: BARNHILL, Stephen, D.; 19 Mad Turkey Crossing. 
Savannah. GA 31411 (US). ZHANG, Zhen; 2055 Middle- 
burg Lane, Mt. Pleasant, SC 29464 (US). 

(74) Agents: TAUTVYDAS, Daiva, K. et al.; Jones & Askew, 37th 
floor, 191 Peachtrec Street, N.E., Atlanta, GA 30303-1769 
(US). 



(81) Designated States: AU, CA, JP, European patent (AT, BE, 
CH. DE, DK, ES. FR, GB, GR. IE, IT, LU, MC. NL. PT, 
SE). 



Published 

With international search report. 



(54) Title: COMPUTER ASSISTED METHODS FOR DIAGNOSING DISEASES 
(57) Abstract 

The present invention is directed to an in vitro, biological fluid-based method for diagnosing, screening or prognosing for diseases in 
which laboratory data are manipulated to obtain a diagnostic index. In one embodiment, a neural network is used to obtain the diagnostic 
index. 



8N500CIO: <WO 9612187A1J 



FOR THE PURPOSES OF INFORMATION ONLY 



Codes used to identify States 
applications under the PCT. 



AT 


Austria 


AU 


Australii 


BB 


Barbados 


BE 


Belgium 


BF 


B tut but Faso 


BG 


.. Bn'.garia 


BJ 


Beam 


BR 


Brazil 


BY 


Belarus 


CA 


Canada 


CF 


Central African Republic 


CC 


Coo jo 


CH 


Switzerland 


a 


Cote d* I voire 


CM 


Cameroon 


CN 


China 


CS 


Czechoslovakia 


CZ 


Czech Republic 


DE 


Germany 


DK 


Denmark 


ES 


Spain 


n 


Finland 


FR 


France 


GA 


Gabon 



party to the PCT on the front pages 



GB 


United Kingdom 


GE 


Georgia 


GN 


Guinea 


GR 


Greece » 


HU 


Hungary 


IE 


Ireland 


IT 


Hah/ 


JP 


Japan 


KE 


Kenya 


KG 


Kyrgyttan 


KP 


Democratic People's Republic 




of Korea 


KR 


Republic of Korea 


KZ 


Kazakhstan 


U 


Liechtenstein 


uc 


Sri Lanka 


LU 


Luxembourg 


LV 


Latvia 


MC 


Monaco 


MD 


Republic of Moldova 


MG 


Madagascar 


ML 


Mali 


MN 


Mongolia 



pamphlets publishing international 



MR 


Mauritania 


MW 


MaUwi 


NE 


Niger 


NL 


Netherlands 


NO 


Norway 


NZ 


New Zealand 


PL 


Poland 


PT 


Portugal 


RO 


Romania 


RU 


Russian Federation 


SD 


Sudan 


SE 


Sweden 


SI 


Slovenia 


SK 


Slovakia 


SN 


Senegal 


TD 


Chad 


TG 


Togo 


TJ 


Tajikistan 


TT 


Trinidad and Tobago 


UA 


Ukraine 


US 


United Slates of America 


UZ 


Uzbekistan 


VN 


Viet Nam 



PCT/US95/01379 



1 

Computer assisted 
methods for diagnosing diseases 

Cross Reference to Related Applications 

The present application is a continuation-in-part of co-pending U.S. 
Application Serial Number 08/315,851 filed on September 30, 1994, which 
is a continuation of U.S. Application Serial Number 07/990,772 filed on 
December 14, 1992, now abandoned, which is a continuation-in-part of 
U.S. Application Serial Number 07/964,486 filed on October 21, 1992, 
now abandoned, which is a continuation of U.S. Application Serial Number 
07/806,980 filed on December 12, 1991, now abandoned. 

Technical Field 

The present invention relates to methods for diagnosing, screening 
or prognosing diseases. More particularly, the present invention relates to a 
method for diagnosing, screening or prognosing diseases in humans or 
animals, and for determining the severity and cause of the disease, by 
measuring blood serum levels of specific, predetermined blood constituents 
and then calculating a diagnostic index based on the relationships between 
those blood constituents. Optionally, demographic information, such as the 
age, race or sex of the patient, as well as the patient's medical history can be 
factored into the diagnosis. 

The present invention further relates to an in vitro, serum-based 
computer assisted method for diagnosing, screening or prognosing 
diseases, utilizing a neural network to obtain a diagnostic index. In 
preferred embodiments of the present invention, the method is used to 
diagnose osteopenia and certain cancers, including but not limited to 
ovarian, breast, testicular, colon and prostate cancer. 
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Background of the Invention 

As used herein, the term "disease" is defined as a deviation from the 
normal structure or function of any pan, organ or system of the body (or 
any combination thereof). A specific disease is manifested by characteristic 
symptoms and signs, including both chemical and physical changes. 
Certain characteristic signs and symptoms can be quantitated by clinical 
chemical analysis to yield important diagnostic information. For purposes 
of this application, the quantifiable signs, symptoms and/or analytes in 
biological fluids characteristic of a particular disease are defined as 
"biomarkers" for the disease. Current diagnostic methods depend on the 
identification and evaluation of these biomarkers, both individually and as 
they relate to one another. 

Generally, the pathological process involves gradual changes that 
become apparent only when overt change has occurred. In many instances, 
pathological changes involve subtle alterations in multiple biomarkers. It is 
uncommon that a single biomarker will be indicative of the presence or 
absence of a disease. It is the pattern of those biomarkers relative to one 
another and relative to a normal reference range, that is indicative of the 
presence of a disease. 

When individual biomarkers do not show a predictable change and 
the patterns and interrelationships among the biomarkers viewed collectively 
are not clear, the accuracy of a physician's diagnosis is significantly 
reduced. Also, as the number of biomarkers significant for a particular 
disease increases, the number of diagnostic patterns increases, 
correspondingly, pattern recognition and accuracy of a physician^ 
diagnosis decreases. >r 

For example, the following biomarkers collectively show 
characteristic changes in the presence of osteoporosis: calcium, phosphate, 
estradiol (follicular, mid-cycle, luteal, or post-menopausal), progesterone 
(follicular, mid-cycle, luteal, mid-luteal, oral contraceptive, or over 60 
years), alkaline phosphatase, percent liver-ALP, and total intestinal- ALP. 
After measuring these biomarkers, a diagnosing physician would next 
compare the measuring to a normal reference range. While some of the 
biomarkers may fall outside the normal reference range, others may fall 
clearly within the normal reference range. In some circumstances, all of the 
biomarker values may fall within a normal reference range. Presented with 
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such data, a physician may suspect that a patient has undergone some bone 
loss, but will be unable to reach a conclusive and meaningful diagnosis as to 
the presence of the disease osteoporosis. 

The characteristic changes in biomarkers associated with some 
5 diseases are well documented; however, the quantitative interpretation of 

each particular biomarker in diagnosing a disease and determining a 
prognosis is not well established. The difficulties inherent in formulating a 
diagnosis from the analysis of a set of laboratory data is best illustrated by 
looking closer at conventional diagnostic methods for a specific disease. A 
10 discussion of the disease osteoporosis follows. 

The term "osteopenia" as used herein means any decrease in bone 
mass below the normal. The term "osteoporosis" as used herein means a 
specific form of generalized osteopenia characterized by a decrease in bone 
density, low bone mass, and microarchitectural deterioration of bone tissue. 
15 Osteopenia encompasses a group of diseases with diverse etiologies 

typified by reduction in bone mass per unit volume to a level below that 
which is necessary for adequate mechanical support. Osteoporosis is the 
result of the gradual depletion of the inorganic portion of the skeleton and 
can be caused by any number of factors. Primary osteoporosis is an age 
20 related disorder that is particularly common in women and is characterized 

by decreased bone mass in the absence of other recognizable causes. 
However, osteoporosis occurs in both men and women. In women it is 
recognized usually at the 5& or 6* decade, following menopause. In men 
osteoporosis is often recognized around their 6* or 7th decade. 
25 The following is a partial list of some of the categories of individuals 

at risk for developing osteoporosis: 

Post-menopausal women 

Cigarette smokers 
30 Heavy users of alcohol 

Users of a variety of drugs, such as steroids 

Female runners and ballet dancers 

Male marathoners consuming two few calories 

Bulimics and anorexics 
35 People with poor diets 

People allergic to dairy products 
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People affected with cancer 

Fair and slim women 

All men and women over the age of 65. 

5 In addition to being female, the three most significant risk factors are 

poor diet, lack of exercise, and being postmenopausal. Other risk factors 
which are associated with osteoporosis include Caucasian or Oriental 
ancestry, a fair complexions, and a family history of osteoporosis. 

The onset of osteoporosis may be insidious or sudden, following 
10 trauma. The most common complaint associated with osteoporosis is back 

pain. Eventually, the pain may spread to the pelvis, the thorax, and the 
shoulders. In the spine, the vertebrae can compress, and the back can take 
on a "bent" appearance. Conditions such as kyphosis (humpback) or 
scoliosis may occur. If the spine becomes deformed, other body parts can 
15 be affected as well. For example, the ribs can be pushed against the pelvis. 

or the stomach can be pushed into the pelvis. In addition to spinal 
problems, osteoporosis can also lead to fractures of the hip, wrist, and nh% 
These fractures can occur with only slight trauma and sometimes with no 
trauma at all. Mazess B., ex al y "Bone Density of the Radius, Spine, anJ 
20 Proximal Femur in Osteoporosis " 7. of Bone and Mineral Research, Vol 

3, pgs. 13-18, (1988); Riggs B. L 9 et a/., "Involutional Osteoporosis", 
New Engl J. Med., Vol. 314 t pgs. 1676-1686, (1986). The changes 
associated with osteoporosis are gradual so osteoporosis is often not 
detected in its early stages. 
25 From the research perspective, osteoporosis can be classified as 

either primary or secondary to another disease. Primary osteoporosis is 
.further classified as juvenile, idiopathic, postmenopausal (Type I) and 
involutional (Type II). It is now understood that accelerated bone loss 
occurs with cessation of menstruation at the time of menopause and in 
30 women who have amenorrhea as a result of prolactin producing pituitary 

tumor, anorexia nervosa, or intense long-distance running associated with 
undernourishment. These situations are all accompanied by estrogen 
deficiency which is likely to be a major determinant of the accelerated bone 
loss. Bone loss also occurs when estrogen therapy is withdrawn. 
35 Bone consists primarily of an extracellular matrix containing (by 

weight) approximately 35% organic and 65% inorganic components. The 
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cells of bone represent a minor component of the bone constituents, yet they 
cany out a major portion of the function of the skeletal system. Bone cells 
help maintain serum calcium concentration within a narrow range to regulate 
mineral homeostasis while also being responsible for the continuous 
5 formation and resorption of the extracellular matrix, allowing response of 

the skeletal system to the mechanical forces resulting from physical activity. 

The major organic component of the extracellular matrix is collagen. 
This protein has a rigid rodlike structure and is composed of three alpha 
10 chains held together in helical fashion by covalent and noncovalent forces. 

Multiple collagen molecules form fibrils, and these fibrils in turn are 
arranged in bundles or fibers. It is the bundles or fibers of collagen that can 
be seen in the light microscope as layers or linear arrays. 

Non-collagen components of bone compose a very small portion of 
15 the organic matrix of the skeleton. This minor fraction consists of proteins, 

glycoproteins, lipids and mucopolysaccharides. Only a few of these 
components have been carefully identified and characterized. Two of the 
proteins that have been isolated and studied are osteonectin and osteocalcin. 
Calcium and phosphorus are the main components of the inorganic 
20 portion of the skeleton. Initially, calcium and phosphorus are deposited as 

amorphous salts but later undergo rearrangements into a crystalline structure 
that resembles hydroxyapatite (Caio(P04>6(OH)2). Several other ions, 
including Na, K, Mg, and CO3 in varying proportions, may be found in the 
skeletal hydroxyapatite. If there has been fluoride intake, there will also be 
25 fluoride in the hydroxyapatite. 

Although considerable effort has been expended in studying the 
mechanism of bone mineralization, it is not fully understood. Several 
theories have been suggested to explain the available data. One theory is 
that calcium and phosphorus ions are present in the extracellular fluid in 
30 amounts exceeding the solubility product of [Ca] x [P]. These ions are 

kept from precipitating by inhibitors of calcification such as pyrophosphate. 
Because osteoblasts contain large amounts of alkaline phosphatase, it has 
been speculated that the activity of this enzyme facilitates mineralization by 
cleaving the phosphate groups, thus altering the Ca:P ratio in the sites of 
35 calcification. 
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In the clinical syndrome of osteoporosis, the reduction in bone mass 
can be attributed to osteopenia because of a dietary deficiency or absorption 
interference of proteins or Vitamin C. It can also be the result of a deficient 
stress stimulus. Osteopenia can also be caused by osteomalacia, a failure of 
proper mineralization of osteoid resulting from bone calcium or phosphorus 
deficiency or both. It can also be caused by insufficient absorption from the 
intestine due to lack of calcium or a resistance to the action of Vitamin D due 
to failure of its conversion to the biologically active forms, 25 
hydroxychole-calciferol and 1.25 dihydroxychole-calciferol, formed by the 
liver and kidney, respectively. In addition, it may well be caused by an 
abnormal rate of osteolysis due to parathyroid hormone stimulation of 
osteoclastic activity (osteocasts are hematopoitic in origin, arising from the 
migration of monocytes to bone). In reality, most cases of osteoporosis, 
when carefully analyzed, reveal evidence of the causes of osteopenia. 

Chemical analysis of blood may reveal calcium, phosphorus, and 
alkaline phosphatase within the normal range. However, an isoenzyme of 
alkaline phosphatase may be significantly increased. Increased bone 
resorption seen in osteoporotic patients, which occurs as a result of the 
action of osteoclasts, usually involves the dissolution of both minerals and 
organic matrix eventually leading to increased excretion of urinary 
hydroxyproline. Serum estradiol which is secreted almost entirely by the 
ovary is significantly decreased in these patients. This observation is further 
corroborated where it has been demonstrated that exogenous estrogen 
therapy in perimenopausal women does delay the onset of post menopausal 
osteopenia. Most experts agree that estrogen appears to decrease bone 
resorption. Weiss N. S. 9 ex aU "Decreased Risk of Fractures of the Hip 
and Lower Forearm with Postmenopausal Use of Estrogen," N. Engl J. 
Med., Vol. 303, pgs. 1195-1198 (1980); Ettinger B., et al % "Long Term 
Estrogen Replacement Therapy Prevents Bone Loss and Fractures," Ann. 
Intern. Med. Vol. 102, pgs. 319-324, (1985). 

It is believed that blood levels of calcium are maintained without 
increasing calcium loss from the bone as long as there are normal amounts 
of estrogen and parathyroid hormone present It is reasonably excepted that 
estrogen antagonizes the effect of parathyroid hormone. Estrogen 
deficiency, as seen in peri- and postmenopausal women, results in an 
increase in the sensitivity of bone to parathyroid hormone. This antagonistic 
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relationship eventually leads to an increase in the resorption of the bone and 
contributes to the development of Type I Osteoporosis. Type II 
Osteoporosis is characterized by reduced calcium absorption, which in turn 
results in increased secretion of parathyroid hormone. Type II Osteoporosis 
occurs in older individuals and is associated with wedge fractures of the 
spine and fractures of the hip. A theory recently advanced by BarnhiU, et al 
suggests that evidence of decreased serum estradiol and increased 
lymphocyte alkaline phosphatase, represents an activated immune system in 
osteopenic postmenopausal women. This finding suggests that "uncovered 
estrogen receptors" may induce an immune reaction which is responsible for 
one form of osteoporosis. Barnhill S.D., et al, "Osteoporosis: A Possible 
Autoimmune Etiology," Ann. of Clin. Lab. Sci. 9 Vol. 17, pgs. 255-256, 
(1987). 

In 1987, Barnhill, etal proposed an autoimmune etiology for some 
forms of osteopenia. In that publication, Barnhill et al showed that 
lymphocyte-derived alkaline phosphatase was present in the blood of 90% 
of severely osteopenic women. See Barnhill, etal, "Osteoporosis: A 
possible Autoimmune Etiology", Ann. of Clin. Lab. Scl, Vol. 17, pgs. 
255-256, (1987)). The concept of a lymphocyte-derived alkaline 
phosphatase is further described by Griffiths, et al See Griffiths, et al, 
"Separation and Identification of Alkaline Phosphatase Isoenzymes and 
Isoforms in Serum of Healthy Persons by Isoelectric Focusing", Clin 
Chem, Vol. 32, pgs. 2171-2177, (1987). Thus, osteoporosis is well suited 
to diagnosis by interpretation of laboratory data. 

Osteoporosis poses a major health problem in the United States, not 
only for those persons who are already affected but for individuals whose 
diet, life style, and body build make them more likely to develop 
osteoporosis as they age. Postmenopausal osteoporosis is a common 
disorder that results in substantial morbidity and mortality. As many as 25 
million American women suffer from osteoporosis. Osteoporosis is 
responsible for approximately 250,000 femoral neck fractures annually in 
the United States, and these fractures are associated with a 20% mortality 
rate within 6 months. The risk is particularly high among the elderly, who 
also tend to lose bone as a result of the aging process. Wu AL, et al, "Bone 
Resorption Rates in Physiological, Senile, and Postmenopausal 
Osteoporosis." / Lab Clin. Med, Vol. 69, pg. 810, (1967). 
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It is estimated from a survey of medical clinics that of those 
individuals living to age 85, 32% of the women and 17% of men will 
fracture hips, weakened by osteoporosis. In addition to the pain and 
suffering caused by these fractures, the monetary cost is great, accounting 
5 for well over 3.8 billion dollars per year in treatment for fractures of 

osteoporosis. Moreover, six to eight months following hip fracture, about 
50% of the osteoporotic patients are in need of assistance with activities of 
daily living and about 25% require nursing home care. Only about 25% of 
the patients fully recover. In view of the costs associated with osteoporosis, 
10 as measured in dollars and in human suffering, osteoporosis has 

increasingly been perceived as a serious and disabling disease, warranting 
substantial involvement on the part of clinical investigators, governmental 
agencies, and pharmaceutical industries to develop and evaluate potential 
treatments and early detection techniques. 
15 The National Institutes of Health 1984 Consensus Conference 

stimulated those interested in bone density with its recommendations for 
"defining persons at risk, and developing safe, effective and low cost 
strategies for fracture protection." See, Office of Medical Applications of 
Research National Institutes of Health, "Osteoporosis: consensus 
20 conference," JAMA, Vol. 252, pgs. 799-802, (1984). Since 1984, there 

has been a dramatic change in the way physicians view osteopenia. With 
the development and popularization of highly sensitive techniques such as 
DPA and quantitative computed tomography, physicians are now capable of 
measuring the density of the proximal femur and the lumbar vertebrae. (For 
25 a general review of the methods currently available to measure bone density, 

see Avioli, LV., ed. f "Metabolic Bone Disease and Clinically Related 
Disorders", W.B. Saunders Company (1990)). 

An early decrease in bone mass can be measured by non-invasive 
assessment of the skeleton by four widely available methods, including 
30 single photon absorptoraetry, dual photon absorptometry, dual-energy x-ray 

absorptometry, and quantitative computed tomography. A brief discussion 
of these diagnostic methods follows. 

A device called a single-photon absorptometer (SPA) is used to 
measure bone mineral content, primarily in the forearm and wrist. The heel 
35 bone can also measured using SPA because the heel bone is thought to be a 

predictor of bone loss in the spine. SPA measures primarily cortical bone, 
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which is also affected by osteoporosis, though not to the same extent as 
trabecular bone. 

The technique of dual-photon absorptoraetry (DPA) provides a 
measurement of the total cortical and trabecular mineral content of the hip 
and spine. DPA uses less radiation than conventional x-rays; but, a scan of 
the spine using DPA still exposes the patient to approximately one tenth of 
the radiation that results from a routine chest x-ray. 

Dual-energy x-ray absoiptometry (DXA) provides a measurement of 
the amount of bone tissue in the hip and spine. This technique is now used 
routinely because it is faster than DPA. 

Unlike DPA and DXA, quantitative computed tomography, more 
commonly called a CAT scan, can measure the density of either bone or just 
the trabecular portion. CAT scans unfortunately expose patients to higher 
doses of radiation than any of the other techniques. 

Radiographic absorptometry (RA) is a method for non-invasive 
measurement of bone mineral x-rays of the hand. Radiographs, taken with 
a standard x-ray machine, are sent to a central laboratory for computer- 
controlled analysis. 

As osteopenia progresses, both the number and thickness of the 
trabecular units decrease, causing fragility of bone and an increased risk of 
fractures. The radiographic manifestations of osteoporosis reflect the 
deficiency of the organic matrix and parallel the gross pathologic findings. 
The cortices are thin and the trabeculae fine and sparse; the skeletal 
structures are therefore more radiolucent than normal. The disease process 
eventually affects almost all of the skeletal structures. The principal areas of 
demineralization arc the spine and pelvis, especially in the femoral neck and 
head. Demineralization is less marked in the skull and extremities. 
Although useful to detect breaks in the bone, ordinary x-rays are not 
sensitive enough to detect osteoporosis until a large amount of bone tissue 
has already been lost, generally from 25% to 40%. By the time 
osteoporosis can be identified by x-ray techniques, the condition is 
advanced. 

Current standard diagnostic techniques, are not effective for early 
detection of osteoporosis. Changes seen in osteoporosis are very gradual, 
and often go undetected in the early stages of the disease. Osteoporosis is 
often not detected in its early stages because bone mass must be decreased 
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by about 30% to 40% before it is apparent using standard x-ray diagnostic 
techniques. Preventing osteoporosis by detecting early bone loss is far 
better than identifying the disease at relatively advanced stages and 
subsequently attempting to prevent its progression. Once major 
5 deterioration has occurred and gaps exist between the ends of fractured 

trabecular, no current treatment can be expected to restore the lost bone. 
Thus, therapeutic efforts must be directed toward prevention and early 
recognition of the progressive disease so treatment can be instituted before 
essentially irreversible structural damage ensues. Cummings S.R., et al t 
iO "Should Perimenopausal Women Be Screened for Osteoporosis?", Ann. 

Int. Med, Vol. 104, pgs. 745-751, (1986); Courpron P., "Bone Tissue 
Mechanisms Underlying Osteoporosis," Orthop. Clin. North Am., Vol. 
12, pg. 513, (1981); Frost H. Af., "Mechanical Determinants of Bone 
Modeling," J.Metabol. Bone. Dis. Rel Res., Vol. 4, pg. 217, (1983). 
15 One of the problems with the current methods for diagnosing 

osteoporosis is that the procedures do not give any information about the 
underlying cause of the osteoporosis, making it difficult to prescribe an 
appropriate course of treatment for the patient. For example, a common 
cause of postmenopausal osteoporosis is an estrogen deficit, which x-ray 
20 techniques cannot measure. Another problem inherent in the current 

diagnostic methods for osteopenia is that all of the current methods require 
expensive, sophisticated medical instrumentation to perform the bone 
density measurements. Additionally, patients must be exposed to x-rays. 
This makes a general screening of high risk populations impractical due to 
25 the expense and unavailability of the necessary instrumentation to the 

average clinic. 

In view of the difficulties associated with extracting a diagnosis from 
the laboratory data for a set of predictive bioraarkers, there is need for 
automated diagnostic systems that are capable of complex pattern 

30 recognition. There have been several attempts at using computational 

models to achieve pattern recognition in diagnostics. One of the most 
popular computational method for making diagnoses from multivariate 
laboratory data has been discriminate function analysis. However, 
diagnostic systems that rely exclusively on classical pattern recognition 

35 technology (geometric, syntactic, template, statistical) are not effective for 

evaluating the characteristic biomarker patterns of many disease states. 
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There is no clear set of rules that accurately describes how to analyze a set of 
biomarkers to reach a diagnosis. 

In recent years, neural networks have been gaining popularity as a 
means for recognizing and analyzing subtle diagnostic patterns in 

5 multivariate laboratory data. Neural networks possess the ability to discern 

patterns and trends too subtle or too complex for humans and conventional 
computational methods to identify. While humans can't easily assimilate 
more than two or three variables at once, neural networks can perceive 
correlations among hundreds of variables. Examples of areas in which 

10 neural networks have been explored for their value in clinical diagnosis 

and/or prognosis include: 

• psychiatry (See Ntulsant, B.H., "A Neural Network as an 
Approach to Clinical Diagnosis", MD Computing, Vol. 7, pp. 25-36 
(1990)); 

15 ♦ autism (See Cohen, L t et al, "Diagnosing Autism: A Neural 

Net-Based Tool", PCAI, pp. 22-25 (May/June 1994); 

• pediatric radiology (See Boone t L.M., et a/., "Neural Networks 
in Radiologic Diagnosis. L Introduction and Illustration", Invest Radiol, 
Vol. 25, pp. 1012-1016, (1990) and Gross, G.W., et al, "Neural 

20 Networks in Radiologic Diagnosis. II. Interpretation of Neonatal Chest 

Radiographs", Invest Radiol, Vol. 25, pp. 1017-1023 (1990)); 

• breast cancer (See Astion, ML., et al, "Application of Neural 
Networks to the Interpretation of Laboratory Data in Cancer Diagnosis", 
Clin. Chem., Vol. 38, No. 1, pp. 34-38 (1992); Yuzheng, W. t et al, 

25 "Artificial Neural Networks in Mammography: Application to Decision 

Making in the Diagnosis of Breast Cancer", Radiology, Vol. 82, pp. 81-87 
(1993); Kappen, H.J., et al, "Neural Network Analysis to Predict 
Treatment Outcome", Annals of Oncology, Vol. 4, Supp. 4, pp. S31-S34 
(1993); and, Ravdin, P.M., et al, " A practical application of neural 

30 network analysis for predicting outcome of individual breast cancer 

patients", Breast Cancer Research and Treatment, Vol. 22, pp. 285-293 
(1992)); 

• ovarian cancer (See Wilding, P., et al, "Application of 
backpropogation neural networks to diagnosis of breast and ovarian 
35 cancer". Cancer Letters, Vol. 77, pp. 145-153 (1994)). 
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• thyroid disease (See, Sharpe, P.K. t et. al\ "Artifical Neural 
Networks in Diagnosis of Thyroid Function from in Vitro Laboratory 
Tests," Clin. Chenu, Vol. 39, No. 11, pps. 2248-2253 (1993)). 

• cervical cancer (See U. S. Patent No. 4,965,725 to Rutenberg)\ 

5 and, 

• cardiology (See U.S. Patent No. 5,280,792 to Leong et al and 
Furlong, J. W. % "Neural Network of Serial Cardiac Enzyme Data: A Clinical 
Application of Artifical Machine Intelligence", Clin. Chenu, Vol. 96, No. 1, 
pp. 134-141 (July 1991). 

L0 Neural networks are capable of pattern recognition particularly suited 

to making diagnoses. Unlike current methods for arriving at a diagnosis 
from a logical set of rules, neural networks do not require explicit encoding 
of process knowledge in a set of rules. Neural networks learn from 
examples. 

15 The present invention provides a new method for diagnosing 

diseases and screening patients at risk for a disease, utilizing laboratory data 
to obtain a diagnostic index. It further provides a method for early detection 
and diagnosis of a disease. As applied to osteoporosis, it generates a new 
level of interest in screening individuals who are at risk for osteoporosis 

20 because this novel diagnostic procedure provides information about the 

underlying cause of the osteopenia. The present invention also provides a 
simple, inexpensive, and rapid method of estimating bone density using 
serum blood levels of specific, predetermined blood constituents. As 
applied to cancer, it provides a simple, inexpensive, accurate and rapid 

25 method for detecting cancer, resulting in earlier and less invasive treatment. 

Summary of the Invention 

The present invention is directed to an in vitro, biological fluid-based 
method for diagnosing, screening or prognosing for diseases in which 

30 laboratory data are manipulated to obtain a diagnostic index. The term 

"biological fluid" includes, but is not limited to, blood, urine, saliva, 
cerebral spinal fluid, synovial fluid, and tears. In the preferred 
embodiments of the present invention, the method of the present invention is 
used to provide a of cancer and osteoporosis. 

35 The present invention is a method for diagnosing, screening or 

prognosing a disease in a human or animal comprising the steps of 
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measuring the concentrations of a predetermined set of biomarkers known to 
be associated with the disease from a biological fluid from the human or 
animal; scaling the digitized values of the analytes; and sending the scaled 
values to a trained neural network, whereby the output values from the 
5 neural network tend toward the upper value when the human or animal has 

the disease and the output values tend toward the lower value when the 
human or animal does not have the disease. 

The present invention also comprises an apparatus for diagnosing, 
screening or prognosing a disease in a human or animal comprising a means 
10 for digitizing the concentrations of a predetermined set of biomarkers known 

to be associated with the disease from a biological fluid from the human or 
animal; a means for scaling the digitized values; and a trained neural network 
coupled to the digitizing and scaling means for generating network output 
values between an upper and lower value, whereby the output values tend 
15 toward the upper value when the human or animal has the disease and the 

output values tend toward the lower value when the human or animal does 
not have the disease. 

The present invention also comprises a method for diagnosing, 
screening or prognosing a disease in a human or animal comprising the 
20 steps of measuring the concentrations of a predetermined set of biomarkers 

known to be associated with the disease from a biological fluid from the 
human or animal, scaling the digitized values of the analytes, and 
introducing the scaled values to a first trained neural network, whereby the 
output values from the first neural network tend toward the upper value 
25 when the human or animal has the disease and the output values tend toward 

the lower value when the human or animal does not have the disease; and 
sending the output value from the first neural network and a second set of 
predetermined biomarkers, which could include one or more of the 
biomarkers in the first set of predetermined biomarkers, to a second trained 
30 neural network, whereby the output values from the second neural network 

tend toward the upper value when the human or animal has the disease and 
the output values from the second neural network tend toward the lower 
value when the human or animal does not have the disease. 

In accordance with the first embodiment of the present invention, a 
35 trained neural network is utilized to determine a diagnostic index 

corresponding to the presence and severity of a disease by analyzing a set of 
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predetermined biomarkers for that disease. In accordance with the 
invention, the concentrations, or in some cases, the presence of certain 
biomarkers related to the incidence of a particular disease are determined for 
a patient These values are then scaled and the scaled values are then sent to 

5 a trained neural network to yield a diagnostic index. A neural network is 

trained by introducing a population of patients in which a disease state is 
known, along with the bioraarker values for those patients and "teaching" 
the neural network to recognize the patterns in the biomarkers. After the 
neural network is trained, biomarker values from patients with unknown 

10 disease states are introduced to the trained neural network. The neural 

network then processes the information to produce a value corresponding to 
a diagnosis of the presence or absence and the severity of a particular 
disease. 

The inventors propose that the artificial neural network, especialh 
15 the multi-layer feedforward network, may, through their weight 

connections, correspond to biomarker patterns that are important for 
categorizing diseases. Additionally, the inventors propose that the neural 
network is able to identify unique patterns of biomarkers associated with a 
variety of disorders that may help to classify borderline cases that do not 
20 appear to fit into either a malignant or benign pattern. 

A second embodiment of the present invention involves a two step 
analysis of the biomarkers by neural network. This avoids the bias created 
by a dominant predictive variable when training a network. The dominant 
biomarker or predictive variable is excluded from the first analysis by neural 
25 network and then included in a second analysis by neural network. For 

example, if age is thought to be the dominant predictive variable in the 
diagnosis of osteoporosis, that variable is not included in the training of the 
first neural network, and the training data set is limited to the other selected 
biomarkers. After obtaining a diagnostic index using the first set of 
30 biomarkers, a second neural network is trained using the diagnostic index 

and the entire set of input variables, including age, to yield another 
diagnostic index. The final diagnostic index is a composition of an artificial 
neural network generated index and results from heuristic analysis using 
other non-numerical patient information. 
35 The present invention further comprises an apparatus capable of 

diagnosing, screening or prognosing a disease comprising a sample 
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receiving means, a sample detecting means capable of detecting the quantity 
of one or more biomarkers in a biological fluid, and, after digitalizing the 
quantity of analyte in the biological fluid, scaling the digitized values and 
introducing the scaled data into a neural network which yields a therapeutic 
5 index on a printer or t optionally, on a video display. 

A third embodiment of the present invention includes a method for 
diagnosing and determining the severity and underlying cause of osteopenia 
using blood concentrations of a specific, predetermined set of blood 
constituents. The method comprises the steps of (1) measuring the severity 
10 of the disease in a set of humans or animals with varying severity of disease 

by a standard method, (2) assigning the severity of disease a numerical 
value on a severity scale, the scale being from no disease to severe disease. 
(3) measuring the blood concentrations of the predetermined set of blood 
constituents in the set of humans or animals with varying severity of the 
15 disease, and then (4) calculating a numerical relationship between the set of 

blood concentrations and the severity of disease. 

The present invention also comprises a simple and rapid method of 
determining severity of osteopenia in a patient In a preferred embodiment 
the method comprises determining the serum level of the following scrum 
20 constituents: calcium, phosphate, total alkaline phosphatase, an alkaline 

phosphatase isoenzyme, estradiol, and progesterone. The alkaline 
phosphatase isoenzyme is preferably Mymphocyte derived alkaline 
phosphatase or blood, liver or intestinal alkaline phosphatase isoenzyme. 
The results of these tests are then introduced into an algorithm. Optionally, 
25 the age of the patient may also be factored into the equation. The bone 

density coefficient that is calculated by the algorithm correlates to a very 
high degree to bone density as measured by standard methods, such as 
radiographic absorptometry, quantitative computed tomography, dual 
photon absorptometry and direct measurement of bone density. The bone 
30 density coefficient that is measured is then compared to an osteopenic 

severity scale. 

Using the six serum constituent concentrations, the present invention 
can be used to determine the osteopenic state of a patient as well as give an 
indication of the underlying cause of the osteopenia. The present invention 
35 can be correlated to any method of measuring bone density simply by 
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recalculating the coefficients in the algorithm using a multiple linear 
regression analysis. 

Accordingly, it is an object of the present invention is to provide a 
method for diagnosing, screening or prognosing and determining the 
5 severity of a disease using the concentrations and values of a predetermined 

set of biomarkers. 

It is a further object of the present invention to provide a more 
accurate method for screening, prognosing and diagnosing breast cancer 
using biological fluid-based laboratory data. 
l0 h is another object of the present invention to provide a more 

accurate method for screening and diagnosing ovarian cancer using 
biological fluid-based laboratory data. 

Yet another object of the present invention is to provide a more 
accurate method for screening, prognosing and diagnosing colon cancer 
15 using biological fluid-based laboratory data. 

Furthermore, it is an object of the present invention to provide a 
more accurate method for screening, prognosing and diagnosing prostate 
cancer using biological fluid-based laboratory data. 

Additionally, it is an object of the present invention to provide a 
20 more accurate method for screening, prognosing and diagnosing testicular 

cancer using biological fluid-based laboratory data. 

Another object of the present invention is to provide a method for 
diagnosing cancer which will provide a better understanding of the probable 
cause of the cancer. 

25 Another object of the present invention is to provide a diagnostic test 

for cancer which can be used to screen large numbers of individuals. 

Another object of the present invention is to provide a simple and 
rapid chemical test for diagnosing osteoporosis. 

Another object of the present invention is to provide a test for 
30 osteoporosis which will also give information as to the underlying cause of 

the osteopenic condition. 

Another object of the present invention is to provide a diagnostic test 
for osteoporosis which can be used to screen large numbers of individuals. 
Yet another object of the present invention is to provide a method for 
35 diagnosing osteoporosis and determining the underlying cause of the 

osteopenia without having to subject the patient to radiation. 
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These and other objects, features and advantages of the present 
invention will become apparent after a review of the following detailed 
description of the disclosed embodiment and the appended claims. 

5 Brief Description of the Figures 

Fig. 1 illustrates a feed forward neural network having multiple 

outputs. 

Fig, 2 illustrates a feed forward neural network having a single 

output. 

[0 Fig. 3 is an equation illustrating the mathematical relationship 

between the input and output of a neuron. 

Fig. 4 is a schematic illustration of the second preferred 
embodiment of the present invention. 

Fig. 5 is a graph showing the change in the neural network 
15 diagnostic index obtained by the method of the present invention over time 

for a ovarian cancer patirenL 

Fig. 6 is a graph showing the change in the neural network 
diagnostic index obtained by the method of the present invention over time 
for a prostate cancer patirenL 
20 Fig. 7 is a graph showing the change in the neural network 

diagnostic index obtained by the method of the present invention over time 
for a colon cancer patirenL 

Fig. 8 is a graph showing the change in the neural network 
diagnostic index obtained by the method of the present invention over time 
25 for a breast cancer patienL 

Fig. 9 shows the patient data used to develop and test the 
osteoporosis neural network diagnostic system. 

Fig. 10 shows the patient data used to test the osteoporosis neural 
network diagnostic system. 
30 Fig. 11 is a program list of the algorithm of the third embodiment 

of the present invention. 

Fig. 12 shows the discrepancies between the DPA measurements 
and the measurements obtained using the method of the present invention 
for individual subjects. 
35 Fig. 13 shows the correlation of the results obtained using the 

present invention and DPA. 
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Fig. 14 shows the normality of values predicted by the present 
invention. 

Fig. 15 shows the normality of residuals (differences between DPA 
measurements and values predicted according to the method of the present 
5 invention), as evidenced by the arrangement of points along a diagonal 

straight line. 

Fig. 16 shows that the variance of residuals. 
Fig. 17 shows the independence of residuals. 
Fig. 18 shows that virtually all observed DPA measurements are 
10 predicted by the method of the present invention. 

Fig. 19 shows the training data used to construct the prostate 
cancer neural network prognostic system. 

Fig. 20 shows the relationships between the input neurons, hidden 
layer and output neuron of a neural network in one embodiment of the 
1 5 present invention. 

Detailed Description 

The first embodiment of the present invention is directed to an in 
vitro, serum-based computer assisted method for screening, prognosing and 
20 diagnosing diseases utilizing a neural network to obtain a conclusive 

diagnosis. The present invention can be adapted to existing diagnostic 
devices that have a collection means, a sample detecting means capable of 
detecting the quantity of an analyte in a biological fluid and a means of either 
printing or displaying the results of the tests on video display means. 
25 The inventors have discovered that biomarkers collectively alter in 

response to a disease process, and collectively constitute a new diagnostic 
biomarker with better disease predictability than the individual biomarkers. 
When the biomarkers are processed and analyzed as a group to yield a single 
diagnostic index, the sensitivity and specificity of the diagnosis is increased, 
30 making it possible for a physician to detect the presence of a disease earlier 

and with greater precision, or estimate a prognosis with greater precision, 
than by analysis of the individual biomarkers. 

In accordance with the method of this invention, blood is first 
collected from a patient, and then centrifuged to separate the red blood cells 
35 from the serum. It should be understood that biological fluids other than 

blood can be used in practicing the present invention. Next, the serum is 
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analyzed, using standard laboratory techniques, to determine the 
concentrations, or in some cases the presence or absence, of specific 
predetermined biomarkers related to the incidence of a particular disease. It 
is to be understood that this process can be carried out automatically in 
5 conventional diagnostic machines. For purposes of illustration, a 

description of the methods for obtaining the values for the biomarkers for 
osteopenia is provided elsewhere in this section. 

The biomarkers relied upon to diagnose a disease by the method of 
the present invention must be predictive of the suspected disease and must 
10 be statistically significant for analysis by a neural network. The selection of 

biomarkers that offers statistically significant discriminating power in the 
diagnosis of disease involves several steps. First an inventory of 
biomarkers that have shown certain relevancy in the diagnosis of the disease 
of interest must be conducted. In general, only the biomarkers that reflect 
15 different aspects of the disease process or other diagnostic information need 

to be included. Second, the selected biomarkers need to have a reasonable 
diagnostic value in terms of sensitivity, specificity, and positive and 
negative predictive powers. The design and implementation of experimental 
protocol from which the biomarkers are developed and evaluated should 
20 also be considered. Third, if the number of candidate biomarkers is large, a 

formal discriminating power analysis may be conducted. However, many 
of the standard statistical analysis methods may not be adequate for highly 
nonlinear classification problems. If the number of candidates are not too 
large, they may be all included in the initial attempt of neural network 
25 training. If one or several of the input biomarkers to the network are 

irrelevant to the classification decision making process, it will be reflected in 
the network connection weights of the trained neural networks. These 
values may then be removed from the biomarker set for a particular disease. 
Other methods for evaluating the statistical significance of a biomarker 
30 selected for analysis by neural network and selecting biomarkers for training 

a neural network are well known in the art 

Biomarkers which meet the criteria delineated above, namely, they 
are predictive of a particular disease and statistically significant for analysis 
by neural network, are identified below for a sample of diseases including 
35 ovarian cancer, prostate cancer, colon cancer, breast cancer, testicular cancer 

and osteoporosis. 
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Ovarian Cancer 


Prostate Cancer 1 


Colon Cancer 


LASA-P® 


I A C A Drift 






PAP 


PA 1Q-Q 


DM/ /UJV 


P^A 


CPA 










Prostate Cancer II 






PAP 






PSA 






CK-BB 






Acid Phosphatase 










Breast Cancer 


Testicular Cancer 


Osteoporosis 


LASA-P® 


LASA-P® 


Calcium 


CEA 


AFP 


Phosphate 


HER2/neu 
in Plasma 


HCG-Beta 


Estradiol 




CA 15-3® 


Progesterone 






ALP 






ALP Isoenzyme 1 






ALP Isoenzyme 2 



5 

A key to the abbreviations used above is provided below: 



AFP: 


Alpha-Fetoprotein 


CA125: 


Cancer Antigen 125 


CA 15-3®** 


Breast Antigens 115D8/DF3 


CA 19-9: 


Carbohydrate Antigen 19-9 


CEA: 


Carcinoembryonic Antigen 


CK-BB: 


Creatinine kinase, BB subfraction 


DM/70K: 


Ovarian marker NB/70K 


HCG-Beta: 


Human Chorionic Gonadotropin, Beta Sub- 
Unit 


HER 2/neu in 
Plasma: 


c-erb B-2 (HER2/neu) oncoprotein in plasma 


LASA-P®*: 


Lipid-Associated Sialic Acid in Plasma 


PAP: 


Prostatic Acid Phosphatase 


PSA: 


1 1 JAM W ^T5SBBi^*nB5il^^^^^^^^^M 



*LASA-P is a registered trademark of DIANON Systems, Inc. 
**CA 15-3 is a registered trademark of Centocor, Inc. 
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A wide number of diseases may be diagnosed in accordance with the 
method of the present invention. To be suitable for diagnosis by the present 
method, biomaikers for the disease must be quantifiable. The biomarkers 
must also be predictive of the disease and must be statistically significant 
5 relative to one another. The method of the present invention is equally 

suited to the diagnosis of any disease in which serum biomarkers can be 
identified, including but not limited to infectious diseases, and genetic 
abnormalities. 

After determining the serum biomarkers for a disease, the biomarker 

10 values are analyzed by a trained neural network to yield a single diagnostic 

value. The most common neural network architecture for pattern 
classification problems is the feedforward network, which typically consists 
of an input layer, one or more hidden layers, and an output layer. Figs. 1 
and 2 illustrate the arrangement of neurons in two different feedforward 

15 networks. 

The elements that make up each layer of a neural network are 
referred to as neurons or nodes. Inputs are fed forward from the input layer 
to the hidden layers and then to the output layer. The number of neurons in 
each layer is determined before the network is trained. Typically, there is 

20 one input neuron or node for each input variable, and one output node for 

each output The inputs to the neural network are predictor variables. 
These predictor variables can be quantitative or qualitative. Neural networks 
make no data distribution assumptions and can simultaneously use both 
quantitative and qualitative inputs. In the present invention, the biomarker 

25 values, after a linear prescaling to values between 0.0 and 0. 1 , constitute the 

input variables. 

The outputs of the network represent output categories. For 
example, a malignancy may be represented by maximal output of the 
malignant output neuron and silence of the benign neuron, whereas a benign 

30 process is represented by maximal output of the benign neuron and silence 

of the malignant neuron. A simple arithmetic function combines outputs of 
the two neurons to yield a single diagnostic index. In the alternative, a 
single output neuron may be used. An output of greater than 0.5 would 
indicate a malignancy and an output of less than 0.5 would indicate a benign 

35 condition. In this way a diagnostic index is directly obtained. 
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The number of hidden layers and the number of nodes in the hidden 
layers are configurable parameters that have a significant influence on the 
performance of the network. In practice, the optimal number of hidden 
neurons is determined empirically. The means for determining the optimum 

5 numbers of hidden neurons is well known to those skilled in the art and 

depends on the complexity of the problem being solved. 

In the present invention, the main network model is a multi-layer 
feedforward perceptron using a backpropogation training algorithm. The 
number of hidden layers and the number of neurons in each hidden layer 

10 was determined to adequately match the level of complexity of the diagnostic 

problem. With the assumption that the samples in the training set are 
representative of all possible situations encountered in real applicationswith 
no significant contradictions, and and the number and stratification of 
samples in the generalization test are statistically adequeate, the criteria 

15 outlined below is used to determine if a chosen network configuration is 

appropriate. 

If the network continues to fail to correctly classify large portions of 
the samples in the training set, even after many adjustments of training 
algorithm parameters, the network complexity should be increased. 
20 On the other hand, if the the network acheives a high rate of 

correctly classifying the training set but fails to accurately classify a large 
number of samples in the testing set, network structure is probably too 
complex for the problem being solved, Le. it has sufficient inherent 
flexibility to fit the training data set, but not sufficient predictive power to 
25 classify the test data set If this is the case, the number of neurons in the 

hidden layers should gradually be reduced, or, if there are mutiple hidden 
layers, the hidden layers should be gradually reduced. 

Note that it is not always necessary to have a large training sample 
set If the samples in a training set have already represented all possible 
30 cases with adequate statistical significance, the addition of new samples 

genewrally does not increase the amount of information in the training 
samples. Instead it may decrease the useful information to nosie ratio in the 
samples. At the other extreme, too small a training set will generally not be 
able to cover all possible variations in the population. The resultant network 
35 often simply memorizes all the cases in the training set and does not 

generalize at all. 
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The input and output layers are not directly connected. Every input 
neuron is connected to every neuron in the following hidden layer and 
neuron in ahidden layer is connected to every neuron in the following 
adjacent hidden layer or output layer, depending on the number of hidden 
5 layers. Each of the multiple connections to a particular neuron is weighted. 

In the hidden and output layers, each node sums the input activations, 
multiplied by the respective connection weights plus a bias term. The 
weighted sum then passes through a non-linear output function, typically a 
sigmoidal function, which gives the network the ability to represent complex 
10 non-linear relationships. A neuron fires if the sum of the weighted inputs to 

it are greater than a threshold value. As illustrated in Fig. 3, once a neuron 
is above a threshold, the magnitude of its output is a sigmoid function of the 
net input The end result of activity in the neural network is the net output, a 
complex nonlinear function of the inputs. 
15 in summary and in accordance with the present invention, first the 

values of the biomarkers for a specific disease are determined and scaled. 
The biomarkers are fed forward from the input layer to the hidden layer (or 
layers) and then to the output layer of the neural network. The number of 
neurons in the input layer is determined before the network is trained and 
20 correspond to the number of biomarkers predictive for a specific disease. 

The biomarkers are preselected. There is one input neuron for each 
diagnostic variable or biomarker, and one output neuron for each desired 
output. Other than the identified biomarkers, diagnostic variables may 
include demographic information such as the age, race or sex of the patient 
25 and whether the patient is pre- or post-menopausal. The number of neurons 

in the output layer depends on the type of output desired. The number of 
neurons in the hidden layer is determined euperically during training. 

The neural network used for diagnosing a specific disease must be 
trained to do so. In accordance with the present invention, the neural 
30 network is trained by back propagation. Back propagation refers to the 

technique of training a neural network to accurately model a set of input and 
output examples by determining the best connection weights between the 
values, and is well known in the art Other techniques which may be used 
to train a neural network for purposes of this invention may include any 
35 other non-linear global optimization technique, such as the genetic search 
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algorithm; however, the feed forward, back propogation network is most 
popular. 

At the beginning of training, the connection weights in the network 
are randomly initialized. The training data is then presented to the network 
one at a time. In accordance with the present invention, the training data 
consists of the biomarker values for a group of patients, and the diagnosis 
for each of those patients. The biomarker values are the input variables used 
to train the network. For each patient, the network uses the patient's 
biomarker values to estimate a diagnosis, which is then compared to the 
actual diagnosis. If the network's diagnosis is correct, then the connection 
strengths and thresholds within the network are not changed, and the next 
patient is presented to the network. If the estimate of the diagnosis is not 
correct, the connection weights and thresholds in both the hidden layer and 
the output layer are adjusted to reduce the size of the classification error. 
After adjustments are made, the next patient is presented. Training proceeds 
until all patients in the training group are correctly classified. If the network 
is not able to completely classify all the data it will train indefinitely, in 
which case training is terminated. 

When training the neural network, the trainer may set the 
decision limits regarding the definition of a classification error, i.e. an 
incorrect diagnosis. The relevant parameter is the error tolerance, which 
specifies how close the estimated output has to be to the actual output to be 
correct. For example, if two output neurons are used and the training 
tolerance is set at 5%, the estimate of malignancy is considered correct if the 
malignant output neuron fires at 95% of maximum and the benign neuron 
fires at 5% of maximum. Similarly, a correct estimate of a benign diagnosis 
means that the benign output neuron firs at 95% of maximum, while the 
malignant neuron fires at 5% of maximum. The methods for determining a 
classification error are well known to those skilled in the art. 

In the preferred embodiment of this invention, if a single output 
neuron is used, a benign normal diagnosis is set at an output of 0.1 and a 
malignant or abnormal diagnosis is set at an output of 0.9. Error tolerance 
is an adjustable parameter and is significant in determining the success of the 
network at making an accurate diagnosis. It is preferred to use either one or 
two output neurons. 
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After the neural network is trained for the desired disease, biomarker 
values from patients with unknown disease conditions are introduced to the 
trained neural network. The neural network then processes the information 
to produce a value corresponding to a diagnosis of the presence or absence 

5 of the particular disease. In accordance with the present invention, this is 

accomplished by using either one single output neuron or two output 
neurons. If two output neurons are used, the output from the two neurons 
are combined to generate a single diagnostic index. 

As illustarted by Fig. 4, in a second embodiment of the present 

10 invention, the diagnostic value obtained by analysis of the biomarkers by a 

trained neural network is further analyzed by a set of heuristic rules in 
combination with additional patient information. The additional patient 
information includes things such as family medical history and demographic 
information. This data is then processed to yield a second single diagnostic 

15 value. 

This invention is further illustrated by the following examples. 

which are not to be construed in any way as imposing limitations upon the 

scope thereof. On the contrary, it is to be clearly understood that resort ma> 

be had to various other embodiments, modifications, and equivalents 
20 thereof, which, after reading the description herein, may suggest themselves 

to those skilled in the art without departing from the spirit of the present 

invention and/or the scope of the appended claims. 

In the following examples which utilize a neural network in the 

analysis of the data, a Neuralshell 2, Release L5 (Ward Systems Group, 
25 Inc.) neural network development program was used for the training of the 

neural network on a Pentium 60 mhz computer (Magitronic, Inc.). 

Example 1 

The following example describes the training of a neural network to 
30 diagnose colon cancer. 

A total of 167 patients ( 48 with diagnosed colon cancer, and 1 19 
normal) were divided into 2 groups, a training set and a generalization 
testing set The training set contained approximately 60% of the patients (90 
patients, 27 with diagnosed colon cance, and 63 normal), and the 
35 generalization testing set contained approximately 40% of the patients (77 

patients, 23 with diagnosed colon cancer, and 54 normal). 
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The initial network architecture was selected based on the level of 
complextity of the classification task. A multi-layer feedforward network 
was used. Selection of the initial architecture involved the selection of the 
number of hidden layers and the number of neurons in each hidden layer. 
5 Several trial iterations were performed to determine and adequeate 

configuration that showed good results on both the training sample set and 
the generalization test sample set This particular neural network had 5 input 
neurons, one hidden layer having 10 neurons and one output neuron. 

Initially, connection weights among the neurons were randomly set 
10 The neural network had five input neurons, corresponding to five input 

variables for diagnosing colon cancer: LASA-P®; CA 19-9, CEA, Age, 
and Sex. During training, the five input variables for each patient were first 
linearly scaled into the continuos range between 0.0 and 1.0. The resultant 
five numbers were then presented as an input vector to the input neurons of 
15 the artificial neural network. 

For each of the input vectors, the network generated an output based 
on the connection weights among the network neurons. The output can be a 
single value or a vector of numbers, depending on the number of output 
neurons used. Each neuron in the network participated in the output 
20 calculation by passing the sum of all inputs to the neuron through a non- 

linear s-shaped function (often a ligistic function) and sending the result to 
each and everyone of the neurons in the following adjacent layer. The 
generated output was compared to the desired "target" output A value of 
0.1 corresponded to a diagnosis of normal and an output of 0.9 
25 corresponded to a diagnosis of of abnormal. The difference was used to 

calculate an error term lo guide the training algorithm, i.e., the 
backpropogation algorithm, in the adjustment of network connection 
weights in an attempt to reduce the differences between network outputs and 
target values over the training sample set. 
30 To counteract over learning, which can occur when a neural network 

is too large for the problem and trains too long, the network was tested on 
the generalization test set every 200 iterations. The networks for which 
generalization testing showed reduced root mean squared error over 
previous tests were saved and used in later calculations or generalizations. 
35 This procedure was repeated until generalization failed to improve over 2000 

successive tests. 
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To accelerate the rate of solution during the initial stage of training, 
the learning rate for the hidden layer was initially set at 0.1 and momentum 
was set at 0.05, for the first several iterations. After the first several 
iterations, the learning rate was lowered to 0.01 and momentum to 0.0005 
to narrow in on the appropriate solution. 

After training, the neural network correctly identified 95% of the 
patients without cancer in the training data set, and 92% of the patients with 
cancer. 

When presented with the generalization test results, the trained 
neural network correctly identified 100% of the people with normal 
diagnosis and 95% of the people with cancer diagnosis. 



15 



20 



25 



Example 2 

The following example illustrates the value of the present invention 
for diagnosing ovarian cancer. 

The biomarker values for LASA-P®, CA 125, and DM/70K were 
measured in a patient suspected of having ovarian cancer. Two of the 
biomarker values fell within the normal reference range, and one value fell 
outside the normal reference range. The values were scaled and analyzed by 
a trained neural network, yielding a diagnostic index of 75, indicating the 
presence of ovarian cancer. The diagnotic index generated by the neural 
network (between 0.1 and 0.9) is multiplied by 100 to obtain a whole 
number. The normal reference range for the diagnostic index was 0-50. 
The biomarker values measured, together with the normal reference range 
for those values are provided below: 



Biomarker 


Patient Value 


Normal Reference 
Range 


LASA-P® 


lfc 


15-20 me/dl 


CA125 


16 


0-35 U/ml 


DM/70K 


45 


0-35 U/ml 


Neural Net 
Diagnostic Index 


75 


0-50 



30 



The patient underwent surgery to treat, the cancer. Fig. 5 graphically 
illustrates the change in the neural network diagnostic index over time, 
before and after surgery. 
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Example 3 

The following example illustrates the value of the present invention 
for diagnosing prostate cancer. 

The biomarker values for LASA-P®, PAP, and PSA were measured 
in a patient suspected of having prostate cancer. All three of the biomarker 
values fell within the normal reference range. The values were scaled and 
analyzed by a trained neural network, yielding a diagnostic index of 90, 
indicating the presence of prostate cancer. The normal reference range for 
the diagnostic index was 0-50. The biomarker values measured, together 
with the normal reference range for those values are provided below 



Biomarker 


Patient Value 


Normal Reference Range 


LASA-P® 


16 


15-20 mg/dl 


PAP 


1.5 


0-2.5 ng/ml 


PSA 


1.0 


0-2.5 ng/ml 


Neural Net 1 90 
Diaenostic Index | 


0-50 



The patient underwent surgery to treat the cancer. Fig. 6 graphically 
illustrates the change in the neural network diagnostic index over time, 
before and after surgery. 



Example 4 

The following example illustrates the value of the present invention 

20 for diagnosing colon cancer. 

The biomarker values for LASA-P®, CA 19-9, and CEA were 
measured in a patient suspected of having colon cancer. All three of the 
biomarker values fell within the normal reference range. The values were 
scaled and analyzed by a trained neural network, yielding a diagnostic index 

25 of 93, indicating the presence of colon cancer. The normal reference range 

for the diagnostic index was 0-50. The biomarker values measured, 
together with the normal reference range for those* values are provided 
below: 
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Biomarker 


Patient Value 


Normal Reference Range 


LASA-P® 


17.9 


15-20 mg/dl 


CA19-9 


5.0 


0-37 U/ral 


CEA 


2.3 


0-2.5 ng/ml 


Neural Net 
Diagnostic Index 


93 


0-50 



The patient underwent surgery to treat the cancer. Fig. 7 graphically 
illustrates the change in the neural network diagnostic index over time, 
before and after surgery. 

5 

Example 5 

The following example illustrates the value of the present invention 
for diagnosing breast cancer. 

The biomarker values for LASA-P®, CEA, HER2/neu in Plasma 

10 and CA 15-3® were measured in a patient suspected of having breast 

cancer. All four of the biomarker values fell within the noimal reference 
range, although two of the biomarker values were on the high end of the 
normal range. The values were scaled and analyzed by a trained neural 
network, yielding a diagnostic index of 91, indicating the presence of breast 

15 cancer. The normal reference range for the diagnostic index was 0-50. The 

biomarker values measured, together with the normal reference range for 
those values are provided below: 



Biomarker 


Patient Value 


Normal Reference Range 


LASA-P® 


18 


15-20 mg/dl 


CEA 


1.0 


0-2.5 ng/ml 


HER2/neu in Plasma 


19 


0-20 U/ml 


CA 15-3® 


8 


0-30 U/ml 


Neural Net 
Diagnostic Index 


91 


0-50 



20 The patient underwent surgery to treat the cancer. Fig. 8 graphically 

illustrates the change in the neural network diagnostic index overtime, 
before and after surgery. 



BNSOOCID: <WO„9612187A1_l. > 



WO 96/12187 



PCTYUS95/01379 



30 

Example 6 

This example illustrates the construction and training of a neural 
network for diagnosis of ovarian cancer. 

5 Table I below provides the data used to train the network (age, 

LASA-P®, CA 125, and DW70K), the actual histological diagnosis, and 
diagnostic index determined by the neural network. A histological index of 
0.1 corresponds to a normal diagnosis and a histological index of 0.9 
corresponds to an abnormal diagnosis. A neural net diagnostic index less 

10 than .5 corresponds to a normal diagnosis. A neural net diagnostic index 

greater than .5 corresponds to an abnormal diagnosis. 

TABLE I 

Training Data for Ovarian Cancer Diagnostic System 



Patient 
Number 




.ASA-P€ 


ZA125 


DM/70K 


II Neural Net 
istological 1 Diagnostic 
Diagnosis fl Index 


94 


34 


15.8 


4 


7 


0.1 | 0.225 


99 


27 


19.1 


32 


7 


0.1 1 0.357 


100 


33 


27 


19 


27 


0.1 1 0.317 


101 


32 


16.8 


7 


0 


0.1 | 0.236 


102 


29 


17.5 


8 


13 


0.1 I 0.186 


106 


25 


18 


3 


0 


0.1 1 0.137 


110 


36 


17.7 


21 


0 


O.l | 0.419 


112 


37 


12.4 


10 


0 


0.1 1 0.318 


114 


27 


16.2 


7 


11 


0.1 1 0.159 


2001 


44 


43.5 


1880 


182 


0.9 | 0.657 


2007 


51 


27 


142 


86 


0.9 


0.75 


2009 


64 


20.4 


55 


32 


0.9 


0.87 


2010 


72 


27.1 


7 


43 


0.9 


0.789 


2015 


39 


22 


60 


55 


0.9 


0.568 


2019 


35 


16.5 


13 


0 


0.9 


0.327 


2020 


66 


23.1 


19 


46 


0.9 


0.79 


2022 


69 


23.2 


189 


28 


0.9 


1 0.905 
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2024 


88 


22.4 - 


47 


32 


0.9 


0.962 


2026 


69 


12.6 


11 




0.9 


0.807 


2027 


72 


25.3 


58 


21 


0.9 


0.923 


2031 


61 


25.1 


1480 


58 


0.9 


0.848 


2032 


79 


24.7 


78 


21 


0.9 


0.946 


2034 


65 


30.6 


439 


45 


0.9 


0.876 


2036 


48 


14.2 


13 




0.9 


0.538 


2037 


72 


16.5 


15 




0.9 


0.858 


2038 


68 


13.8 


25 


11 


0.9 


0.851 


2039 


69 


12.1 


10 


3 


0.9 


0.796 


2040 


79 


25.2 


50 


45 


0.9 


0.939 


46 


42 


22.6 


3 


15 


0.1 


0.34 


49 


25 


16.1 


6 


7 


0.1 


0.14 


51 


28 


15.2 


17 


0 


0.1 


0.251 


52 


40 


19.1 


0 


0 


0.1 


0.303 


53 


27 


25 


4 


5 


0.1 


0.17 


56 


31 


24.6 


0 


0 


0.1 


0.196 


57 


33 


13.7 


7 


3 


0.1 


0.235 


58 


25 


24.8 


12 


9 


0.1 


0.19 


59 


25 


25 


9 


11 


0.1 


0.169 


60 


25 


17.4 


5 


0 


0.1 


0.146 


61 


35 


23.3 


0 


15 


0.1 


0.217 


62 


29 


15.4 


12 


0 


0.1 


0.229 


63 


28 


18.3 


2 


1 


0.1 1 0.16 


65 


33 


15.4 


8 


9 


0.1 R 0.236 


68 


30 


22.8 


2 


15 


0.1 


0.171 


70 


28 


15.4 


24 


1 


0.1 


0.305 


71 


31 


20 


6 


19 


0.1 


0.193 


73 


30 


18 


9 


0 


0.1 


0.228 


75 


35 


23.1 


1 


0 


0.1 


0.25 


76 


35 


20.4 


3 


0 


0.1 


0.258 


80 


34 


18.9 


1 


6 


0.1 


0.215 


82 


32 


18.2 


72 


152 


0.1 


0.435 
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84 


33 


17.4 


19 


21 


0.1 


0.305 


85 


26 


15.6 


7 


5 


0.1 


0.156 


89 


26 


20.7 


12 


0 


0.1 


0.205 


92 




zu. / 


1 


u 


0.1 


0.164 


2042 


68 


16.4 


22 


3 


0.9 


0.852 


2043 


63 


14.5 


9 


0 


0.9 


0.733 


2044 


84 


17.6 


23 


17 


0.9 


0.931 


2046 


72 


17.8 


8 


15 


0.9 


0.813 


2047 


61 


13.4 


146 


15 


0.9 


| 0.852 


2048 


44 


20.8 


21 


0 


0.9 | 0.567 


2049 


60 


17.8 


29 


0 


0.9 1 0.822 



An evaluation of the training data set (n=61) yielded a sensitivity of 96<7< 
and a specificity of 100%. 

Table II provides the data used to test the network trained with the 
5 data provided in Table I. 



TABLE II 

Testing Data for Ovarian Cancer Diagnostic System 



*alient 
Number 


Age 


.ASA-P® 


CA 125 


DW70K 


Histological 
Diagnosis 


Neural Nei 
Diagnostic 
Index 


103 


43 


14.1 


8 


15 


0.1 


0.371 


104 


28 


15.9 


11 


0 


0.1 


0.211 


105 


46 


15.9 


15 


0 ' 


0.1 


0.529 


107 


44 


23.6 


10 • 


26 


0.1 


0.416 


108 


32 


21.1 


9 


3 


0.1 


0.258 


109 


42 


14.1 


13 


12 


0.1 


0.407 


111 


35 


19.3 


9 


6 


0.1 


0.29 


113 


39 


18.7 


4 


7 


0.1 


0.305 


115 


31 


21.8 


10 


0 


0.1 


0.259 


2002 


69 


16.4 


29 


21 


0.9 


0.869 


2003 


53 


13 


11 


0 


0.9 


0.599 


2004 


58 


34.3 


80 


90 


0.9 


0.823 
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9<VK 
L\J\JD 


64 


17 R 
i / .o 


76 


12 


0 Q 


u.ooz 


OAAA 
ZVA/O 


** / 






104 




A AQQ 


OAAfi 
ZUUO 


^A 


1Q 


A 
D 




A Q 


A SAQ 


0A1 1 
ZU1 1 


83 


14 i 

14. J 


1 A 
ID 




A o 

v.y 


A o i o 

U.V1 / 


oai O 
ZU1Z 


AA 


04 1 
Z4..3 


OA 
ZU 


JO 


A o 


A OOO 
U. 1 55 


oai i 

ZU13 


CI 
ol 


01 0 

Zl . / 


10 
1Z 


1 A A 


A o 


A fifl 
U.oo 


OAI A 
ZU14 


AO 




1 0 
1Z 


3V 


A O 


A OQC 


2016 


55 


17.3 


>1A 

4y 


5 


A A 

0.9 


A 1 

0.521 


2017 


68 


19.6 


11 


0 


A A 

0.9 


A O 1 o 

0.813 


2018 


56 


17.8 


V4 


3 


A A 
0,V 


A OOO 

0.83/ 


oao 1 

2021 


68 


12.2 


5 


0 


A A 
0.9 


A T^C 1 

0.761 


2023 


A C 

45 


1 C A 

15.9 


71 


1 


A A 
0.9 


0.713 


2025 


61 


t *T A 

17.9 


34 


5 


A A 

0.9 


A O A a 

0.846 


2028 


43 


17.6 


30 


0 


0.9 


0.617 


2029 


42 


16 


16 


0 


0.9 


0.47 


2030 


51 


21 n 


820 


57 


A A 

0.9 


A *7C 

0.75 


47 


31 


21.7 


5 


a 
0 


A 1 
0. 1 


A OOO 

0.ZZ3 


vIO 

45 


5y 


1 O 1 

17.1 


3 


1 


A 1 

u. 1 


A 1(\A 


CA 

50 


OA 

30 


1 1 O 

17.2 


5 


OA 

Zo 


A 1 
0. 1 


A 1 A 
0. lO 


C A 

54 


44 


O 1 1 

21.1 


>• 
4 


0 


A 1 
0. 1 


A A 1 1 
U.41 1 


55 


OA 
ZD 


IB. J 


1 
1 


*> 

z 


A 1 
U. 1 


U. 1 5D 


AvI 
<>4 


O^ 
ZD 


1*7 0 
1 /.Z 


o 
z 


A 


A 1 


A 1 7 1 


AA 
DO 


45 




c 

J 


1 
1 


A 1 




AO 
O/ 


00 
Z / 


1 A 7 
If./ 


o 


A 
U 


0 1 

v. 1 


0 14Q 




29 


17 2 


3 


12 


0.1 


0.158 


72 


26 


20.9 


2 


4 


0.1 


0.142 


74 


31 


19.8 


3 


20 


0.1 


0.173 


77 


25 


15 


10 


0 


0.1 


0.168 


78 


25 


12.4 


4 


0 


0.1 


0.132 


79 


35 


20.7 


1 


6 


0.1 


0.233 


81 


26 


17.2 


9 


0 


0.1 


0.178 


83 


45 


17.1 


2 


10 


0.1 


0.372 


86 


32 


15.6 


11 


0 


0.1 


0.262 
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87 


31 


21.5 


6 


11 


0.1 


0.21 


88 


26 


17.2 


3 


0 


0.1 


0.145 


90 


37 


16.5 


2 


2 


0.1 


0.262 


91 


30 


16.9 


7 


0 


0.1 


0.211 


93 


34 


21.7 


8 


1 


0.1 


0.285 


95 


35 


20 


13 


0 


0.1 


0.339 


96 


34 


22A 


7 


i 
J 


U. 1 




97 


37 


15.9 


6 


6 


0.1 


0.284 


98 


38 


23.5 


38 


0 


0.1 


0.624 


2033 


72 


17.6 


3 


0 


0.9 


0.806 


2035 


38 


21.6 


19 


12 


0.9 


0.422 


2041 


59 


21.5 


1660 


9 


0.9 


0.857 


2045 


61 


23 


436 


46 


0.9 


0.844 


2050 


55 


28.4 


480 


139 


0.9 


0.794 



The testing data (n=59) yielded a sensitivity of 92% and specificity 
of 91%. 

The combined results of the training and testing data yielded a 
5 sensitivity of 94% and specificity of 96%. 

Fig. 20 provides the equations describing the weight matrices 
between the layers in the neural network 

Table III below compares the sensitivity and specificity of a 
diagnosis obtained by the method of the present invention with the 
10 sensitivity and specificity of each individual biomarker. As is evident from 

the results, the individual biomaikers may be combined by this method to 
yield a 'single diagnostic index having greater sensitivity than the individual 
tests for purposes of diagnosis. 
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TABLE IU 

Comparision of Ovarian Cancer Diagnostic Methods 



BioMarker Test 


Sensitivity 


Specificity 


Carbohydrate Antigen 125 
(CA 125) 


74%* 


99%* 


Dianon Marker 70K 
(DW70K) 


68%* 


99%* 


Lipid Associated Sialic Acid 
(LASA-P®) 


79%* 


96%* 


Neural Network Diagnostic 
Index for Ovarian Cancer 


94%** 


96%** 



5 * BioStatistics provided by Dianon Systems, Inc. 

** BioStatistics based on initial pilot study performed by Horns 
Therapeutics, Inc. 

Example 7 

10 

This example illustrates the construction and training of a neural 
network for diagnosis of osteoporosis. 

Fig. 9 provides the data used to develop, i.e. train and test, the 
neural network to diagnose osteoporosis. The biomarkers selected included 
15 age, calcium, phosphate, estradial (ETWO), progesterone, total alkaline 

phosphatase, total intestinal alkaline phosphatase, and % liver alkaline 
phosphatase. Fig. 9 further includes the diagnostic index obtained by the 
neural network. 

Fig. 10 provides the data use to test the network trained with the data 
20 in Fig. 10, and the neural network diagnostic index obtained. 

The second embodiment of the present invention also relates to a 
method for diagnosing diseases in humans or animals using blood 
concentrations of specific, predetermined blood constituents. More 
specifically, the present invention relates to a method for detecting the 
25 occurrence of osteopenia, diagnosing osteoporosis and determining its 

severity and underlying cause. The present invention also facilitates the 
periodic monitoring of specific physiological functions which may indicate 



BNSOOCID: <WO 96121B7A1 J_> 



WO 96/12187 



PCTYUS95/01379 



36 



the onset of osteopenia and correlates to bone mineral density measurements 
determined by various standard methods. 

In practicing one aspect of the present invention, the severity of 
disease in a set of humans or animals with varying severity of disease is 
5 measured by a standard method or methods. The measurement is then 

assigned a numerical value corresponding to a severity scale. The scale 
ranges from humans or animals with no disease, to humans or animals with 
severe disease. The scale is preferably a numerical scale. For example, one 
could assign a value which corresponds to normal or slight disease, another 
10 value which corresponds to moderate disease and a third value which 

corresponds to severe disease. 

The concentration of a predetermined set of blood constituents in the 
set of humans or animals with varying severity of disease is then 
determined. According to the present invention, it is preferable to measure 
15 the blood constituents in the same set of humans or animals in which the 

severity of disease was measured by the conventional method or methods. 

After determining the blood concentrations of the predetermined set 
of blood constituents in the set of humans or animals with varying severity 
of the disease, mathematical manipulations are performed in which the 
20 numerical value obtained using the conventional diagnosis is set to equal the 

result of a mathematical model utilizing the set of blood constituent 
measurements. For example, the relationship can be made using a multiple 
linear regression analysis. It is to be understood that other mathematical 
models could be used to determine a correlation between the serum 
25 constituent concentrations and a standard method of determining severity of 

the disease. The concept of using mathematical nndels to determine the 
relationship or correlation is considered to be part of the present invention. 
Standard statistical analyses that are well-known to those of ordinary skill in 
the art are then performed to determined the confidence level of the 
30 correlation between the diagnosis by conventional means and the set of 

blood constituents. These statistical analyses include chi-square tests. 

An example of practicing one embodiment of the present invention is 
a method for diagnosing osteopenia in a human or animal. The method 
preferably utilizes six blood constituents. These constituents are calcium, 
35 phosphate, total alkaline phosphatase, an alkaline phosphatase isoenzyme, 

estradiol, and progesterone. The alkaline phosphatase isoenzymes preferred 
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for practicing the present invention include lymphocyte-derived alkaline 
phosphatase isoenzyme and bone, liver or intestinal alkaline phosphatase 
isoenzymes. The present invention includes calculating a bone density 
quotient using the aforementioned six blood constituents by entering the 
5 values for the tests into an algorithm that is calculated using a multiple linear 

regression analysis. Optionally, the age of the patient may also be 
incorporated into the equation. The bone density coefficient derived from 
the algorithm allows one to diagnose the osteopenic state of the patient, 
including the severity of the disease. 
10 In addition to diagnosing the osteopenic state of the human or 

animal, an indication of the underlying cause of the osteopenia can be 
determined using the present invention. For example, by practicing the 
present invention as described herein, one can determine whether the 
osteopenia in a human or animal is caused by post-menopausal lack of 
15 estrogen or is caused by some other condition, such as cancer. This allows 

the attending physician to be better able to prescribe the appropriate 
treatment for the osteopenia. 

Five of the serum tests that are used in the present invention are tests 
that are commonly performed by clinical laboratories. The test for t- 
20 lymphocyte derived alkaline phosphatase is experimental only; however, the 

test for blood, liver and intestinal alkaline phosphatase isoenzymes are also 
known. The type of test used to determine the six serum constituents is not 
critical to the present invention as long as the tests give accurate blood 
concentrations of the constituents being measured. 



25 



Calcium & Phosphorus: 

Maintenance of calcium and phosphorus homeostasis involves the 
participation of 3 major organs: the small intestine, the kidney, and the 
skeleton, and is regulated by various hormones. Calcium enters the body 

30 through the diet and is absorbed into the circulation from the small intestine. 

Calcium absorption occurs by two processes, active transport and passive 
transport. Approximately 98% of the calcium and 85% of the phosphorus 
in the adult is present in the skeleton primarily as hydroxyapatite, which is a 
crystal lattice compound of calcium, phosphorus, and hydroxide. 

35 Remaining calcium is present in extra cellular fluid, some types of tissue, 

and skeletal muscle. Phosphorus is combined with lipids, proteins, 
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carbohydrates, and other organic substances. Of critical importance to 
calcium homeostasis is the fact that less than 1% of total skeletal reservoir of 
calcium is rapidly exchangeable with extracellular fluid. In addition to its 
obvious importance in skeletal mineralization, it is also vital for blood 
coagulation, neuromolecular condition, maintenance of normal tone and 
excitability of skeletal and cardiac muscle, and preservation of all membrane 
integrity and permeability, particularly in terms of sodium and potassium 
exchange. 

Free or ionized calcium accounts for 50% of total calcium. About 
5% of total calcium is complexed with a variety of anions, particularly 
phosphate and citrate. The remaining 45% of calcium is bound to plasma 
proteins. Both ionized calcium and calcium complexes are freely dialyzable. 
Acidotic and alkalotic condhions adversely affect the ionized calcium level in 
the blood. In metabolic bone disease such as hyper- or 
hypoparathyroidism, Pagets of bone, Vitamin D deficiency, and renal 
osteodystrophy, calcium and/or phosphorus levels are significantly altered. 

Serum Calcium Determination: 

In one method of practicing the present invention, the procedure for 
the determination of calcium is based on the interaction of the chromogenic 
agent o-cresolphthalein complexone (Sigma Diagnostics Calcium Agent, 
Sigma Chemical Co., St Louis, MO ) which complexes with calcium cation 
in an alkaline medium to form a purple colored complex which has an 
absorbance maximum at 575 nm. The intensity of the color measured at 575 
nm is directly proportional to the calcium concentration in the given sample. 
The o-cresolphthalein complexone contains 8-hydronyquinoline which 
prevents interference from magnesium ions. Robertson, W.G., et al., 
"Calcium Measurements In Serum and Plasma - Total and Inonized," Crit. 
Rev. Clin. Lab. Sri.. Vol. 11. pg. 271, (1979); Sigma Diagnostics Product 
Insert: Calcium, Procedure No. 587, Reissued May 1989. The calcium 
level in serum is in the range of about 9.2 to 1 1.0 rag/dl which is composed 
of 3 distinct fractions. 
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Specimen Collection: 

In one method of practicing the present invention, serum or 
heparinized plasma are suitable for the test. Anticoagulants other than 
heparin should not be used. 

Procedure: 

Sigma Diagnostics Calcium Reagent is used with Abbott Spectrum 
High Performance Diagnostic System for the quantitation of calcium in 
serum. It is generally considered good laboratory practice to run a 
calibration, linearity, and a quality control prior to the determination of 
patient samples. Calcium levels up to 16 mgs % can be measured by this 
method. 

Serum Phosphorus Determination: 

In one method of practicing the present invention, the procedure for 
the determination of serum phosphorus is based on the interaction of 
inorganic phosphorus with ammonium molybdate in the presence of sulfuric 
acid. Barnhill S„ et aL % "Osteoporosis: A Possible Autoimmune Etiology," 
Ann. of Clin. Lab. ScL. Vol. 17, pg. 255, (1987). This reaction produces 
unreduced phosphomolybdate complex. The absorbance of this complex at 
340 nm is directly proportional to the inorganic phosphorus present in the 
given sample. (Sigma Diagnostics Product Insert: Phosphorus, Inorganic, 
Procedure No. 360-UV. Previous Revision December 1985. Reissued 
July (1988)). Most of the phosphorus in extracellular fluid is inorganic and 
in the range of 2.4 to 4.7 mg/dl. 

Specimen Collection: 

Serum or heparinized plasma are preferred for the test. 
Anticoagulants other than heparin should not be used. 
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Procedure: 

Sigma Chemical Company's phosphorus inorganic reagent can be 
used with the Abbott Spectrum High Performance Diagnostic System for 
qualification of the inorganic phosphorus in serum. It is customary to 
5 perform calibration, linearity, and quality control studies prior to the 

determination of test samples. Phosphorus levels up to 12 mgs % can be 
measured by this procedure. 

Total Alkaline Phosphatase: 
10 The demonstration that bone is rich in alkaline phosphatase (ALP) 

and that normal serum contains the same or a similar enzyme led to the study 
of serum ALP levels in patients with disease of bone. See, Courpron P., 
"Bone Tissue Mechanisms Underlying Osteoporosis," Orthop. Clin. North 
Am., Vol. 12, page 513, (1981). This is especially true in osteitis 
15 deformans, hyperparathyroidism and bone neoplasm. ALP is also elevated 

in hepatic disease; however, this can be distinguished by other corroborative 
laboratory procedures and clinical features. In some situations, as in 
osteoporosis, ALP per se may not be above the reference range; however, 
the isoenzyme of ALP is increased. It is well understood that total serum 
20 ALP in normal subjects consists of isoenzymes contributed from liver, 

bone, renal, pulmonary, placental and intestinal sources, among others. 
The isoenzymes and isoforms of various tissue origins can be further t 
separated and visualized by an Isoelectric Focusing Technique. Estimation 
of total serum ALP by isoelectric focusing techniques offers great potential 
25 in the investigation of metabolic bone diseases. 

Serum Alkaline Phosphatase Determination: 

In one method of practicing the present invention, serum ALP 
activity can be measured using various phosphate esters as substrates. 
30 Sigma Chemical Company's alkaline phosphatase reagent measures serum 

ALP activity by a kinetic method. The reagent for the test contains p- 
nitrophenyl phosphate, carbonate buffer, magnesium ions and mannitol. 
Mannitol present in the reagent acts as a phosphate acceptor during the 
enzyme reaction. McComb R.B., et a/., Alkaline Phosphatase. Plenum 

35 Press, New York, (1979); Gundbery, Af., Alkalinr Phosphatase and 

Osteocalcin, Primer on M etabolic Bone Phases and Disorders of Mineral 
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Metabolism 1st Edition, Published by Am. Soc. For Bone and Mineral Res. 
Editor-Murry J, Favus, pgs. 74-76, (1990); Sigma Diagnostics Product 
Insert: Alkaline Phosphatase (ALP), Procedure No. 245, Reissued March 
1989; Kaplan TA.,et al, Clinical Chemistry. St Louis, e.v. Mosby 
5 Company, (1987). 



Specimen Collection: 

Preferably, serum or heparinized plasma is used for the total alkaline 
phosphatase and alkaline phosphatase isoenzyme determination tests. 
10 EDTA, oxalate, citrate and fluoride are inhibitors of ALP and are not 

suitable anticoagulants. 



Procedure: 

Alkaline phosphatase (ALP) diagnostic reagent manufactured by 
15 Sigma Chemical Company is used with the Abbott Spectrum High 

Performance Diagnostic System for the quantitation of inorganic 
phosphorus in the serum. Serum ALP hydrolyzes p-nitrophenyl phosphate 
to p-nitrophenol and inorganic phosphate. The hydrolysis occurs at alkaline 
pH and the absorbance is directly proportional to the ALP activity of the 
20 serum sample. As is the usual practice, calibration, linearity and quality 

control assays should be performed prior to the determination of test 
samples. 

ALP levels up to 1200 U/L can be measured by this procedure. 
Normal ranges are: 
25 Infants 50-165 (U/L) 

Adult 20-70 (U/L) 

Child 20-150 (U/L) 

60 years 30-75 (U/L) 



30 Alkaline Phosphatase Isoenzyme: 

Alkaline phosphatase isoenzyme levels have been primarily used to 
aid in the differential diagnosis of liver and bone disorders and also to 
indicate placental growth patterns. Isoelectric focusing of alkaline 
phosphatase isoenzyme on a precast gel gives 12 discernible bands of 
35 different organ and tissue origin. Of these 12 Bands, Band 10, appearing at 

4.73 isoelectric point, is of placental origin. It is released by placental 
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syncytiotrophoblast and/or macrophage, and is T-lymphocyte mediated. In 
the alternative, samples for ALP isoenzyme can be analyzed by 
electrophoresis methods which identify mainly liver, bone, and intestinal 
isoenzymes. The Ciba Corning Alkaline Phosphatase Isoenzyme (ALP) 
System may be used for the qualitative and quantitative determination of 
ALP isoenzymes in human serum by electrophoresis. 

The liver isoenzyme is the most frequently encountered alkaline 
phosphatase isoenzyme in pathologically increased serum alkaline 
phosphatase activity. Increased liver isoenzyme is encountered in a variety 
of liver and hepatobiliary diseases. These include hepatic cirrhosis, primary 
biliary cirrhosis, congestive cirrhosis, and hepatic carcinoma. Elevation of 
the liver isoenzyme is also a sensitive indicator of cholestasis and liver 
infiltration. 

Another of the alkaline phosphatase isoenzymes that can be used in 
practicing the present invention is lymphocyte-derived alkaline phosphatase 
The preferred method of measuring this isoenzyme is by isoelectric focusing 
electrophoresis. A commercially available isoelectric focusing apparatus that 
is capable of measuring alkaline phosphatase isoenzymes is made by 
IsoLab. Inc. (Resolve® -ALP. IsoLab, Inc., Akron, Ohio). On the 
Resolve®-ALP isoelectric focusing apparatus, the alkaline phosphatase 
isoenzyme that is used in the present invention resolves as band 10 in the 
electrophoretic pattern. In the preferred bone density algorithm defined 
hereinbelow, the value for the band 10 alkaline phosphatase isoenzyme was 
assigned 0 if the band were missing or very weak, and 1 if the band were 
present The electrophoresis gels can also be scanned with a densitometer 
and more quantitative values can be assigned to the alkaline phosphatase 
isoenzyme concentration. As used herein, the band 10 alkaline phosphatase 
isoenzyme shall be designated lymphocyte-derived alkaline phosphatase 
isoenzyme, or I-Alkp. 

It is to be understood that other methods of measuring I-Alkp could 
be used. These methods include, but are not limited to, enzyme-linked 
immunoassay techniques (ELISA), radioimmunoassay techniques, affinity 
columns, and isoelectric focusing columns. It is important to note that 
measurement of band 10 on the Resolve®-ALP apparatus does not correlate 
at all with osteopenia in general. There are many other abnormal conditions 
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which result in a higher than normal band 10 lymphocyte-derived alkaline 
phosphatase isoenzyme (see Barnhill, et cU. and Griffith, et a/., supra.). 

A number of pathological conditions can lead to an elevated level of 
bone alkaline phosphatase isoenzyme. The highest levels of bone 

5 isoenzymes are usually found in Pagets disease. Increased levels are also 

encountered in rickets, bone cancer, osteomacacia and celiac sprue. Renal 
disorders can also result in increased levels. These include renal failure, 
primary hyperparathyroidism, secondary hyperparathyroidism induced by 
long term hemodialysis and malabsorption. 

10 Increased levels of intestinal alkaline phosphatase are encountered in 

a variety of diseases of the digestive tract These include intestinal infection 
and ulcerative lesions of the stomach, duodenum, small intestine and colon. 

Procedure: 

15 in the preferred method of practicing the present invention, the 

alkaline phosphatase isoenzymes are separated by electrophoresis in a 
buffered agarose system. After electrophoresis the isoenzymes are detected 
by incubating the gel with a fluorescent compound such as 4- 
methylumbelliferyl phosphate. The fluorescence formed during the 

20 interaction is quantitated using a Ciba-corning 7 10 densitometer at 385 nm 

Total alkaline phosphatase activity of the patient is required to interpret the 
liver, bone and intestinal ALP isoenzymes. 

Estrogen and Progesterone: 

25 The estrogens are steroids that have a ring containing three 

unsaturated double bonds. The ovary, as well as the testes and adrenal 
gland, has the capacity to synthesize estrogens from androgens, 
androstenedione and testosterone. During the follicular phase of the 
menstrual cycle, ovarian secretion represents only one third of total estrogen 

30 production. In contrast to estradiol, which is secreted almost entirely by the 

ovary, most estrone is derived from peripheral conversion of 
androstenedione and from estradiol metabolism. During menopause, 
estradiol concentrations steadily decrease to approximately 15 percent of 
premenopausal levels. In healthy postmenopausal women, the ovaries do 
35 not secrete significant quantities of estrogens, and virtually all estrogen 



BNSDOCID: <WO 9612187A1.I > 



PCTAJS95/01379 



44 

produced is from peripheral conversion of androstenedione made by the 
adrenal. Estradiol is present in the serum as follows: 



Follicular Phase 




Early 


30-100 ng/L 


Late 


100-400 ng/L 


Luteal Phase 


50-150 ng/L 


Postmenapausal 


<20 ng/L 



In menstruating females, progesterone is secreted mainly by the 
corpus luteum of the ovary. It is partially responsible for cyclic changes in 
the endometrium that are necessary for attachment and growth of an 
embryo. Progesterone levels are low prior to the mid-cycle gonadotropin 
surge. Shortly after the gonadotropin surge, they begin to rise rapidly, 
reaching peak levels during the middle of the luteal phase. Thereafter, a 
progressive fall occurs with barely detectable progesterone levels prior to 
menses. Although progesterone in large amounts produces a negative 
feedback on gonadotropin secretion, it is not the major component in the 
negative feedback system of ovarian steroids. Function of the corpus 
luteum can be assessed by measuring serum progesterone concentration. 
Progesterone is present in the serum as follows: 



Follicular Phase 


0.1 - 1.5 ng/L 


Luteal Phase 


2.5 - 28.1 ng/L 


Mid-Luteal Phase 


5.7 - 28.1 ng 


Over 60 Years 


0.0 - 0.2 ng/L 



Procedure: 

In a preferred method of' practicing the present invention, progesterone 
measurements can be obtained by a radioiramuno-assay method using 
antibody coated tubes (Diagnostic Products, Los Angles, California). 
Estradiol measurements can be performed by a Microparticle Enzyme 
Immunoassay (IMX), available from Abbott Diagnostics, Abbott Park, IL. 



WO 96/12187 



PCT/US95/01379 



45 



Bone Mineral Density (BMD) Measurements 

The conventional methods of diagnosing osteopenia which may be 
used when practicing the present invention are outlined below: 

Dual-Photon absorptometry (DP A) is widely used to assess bone 
5 mineral content and bone mineral density. Johnston Conrad C, et ai f 

"Clinical Use of Bone Densitometry", New Eng. J. of Med., Vol. 324, pgs. 
1 105-1 109, (1991). DPA uses transmission scanning with photons from a 
radioisotope source, such as 153 Gd, that emits two energy peaks, thus 
allowing bone density to be measured independent of soft tissue. DPA 
10 measurements are performed on lumbar spine (L1-L4) and femur (femoral 

neck, Ward's Triangle and trochanter) and the average determined 
separately. The overall average of both hip and spine can also be 
determined. BMD can be measured by DPA using - Gadolinium as the 
source (Lunar DP3, by Lunar Radiation Corporation, Madison, Wisconsin). 
15 The present method employs transmission scanning using 44 and 

100 KeV photon energies from a one Ci 153 Gd source to allow computation 
of the mineral content of bone independent of soft tissue thickness. Bone 
mineral density, expressed in g/cm2, is derived by dividing bone mineral 
content (BMC) by the projected area of the scanned bone. Peppier W.W., 
20 et a/., "Total Body Bone Mineral and Lean Body Mass by Dual-Photon 

Absorptometry: L Theory and Measurement Procedure, * CALCIP, Issue 
Int.. Vol. 33, pg. 353, (1981); Shipp G. et aL, "Precision of Dual-Photon 
Absorptometry," CALCIP, Issue lnu f Vol. 42, pgs. 287-292, (1988). 
During spine scans, the detector moves in a rectilinear pattern at a rate of 
25 5mm/sec and with scan lines 4.5 mm apart. The BMC and BMD are 

calculated in lumbar vertebrae 1 through 4 (including intervertebral discs) 
with a software version supplied by Lunar Radiation Corporation of 
Madison, Wisconsin. 

During femur measurements, the scanner moves at a rate of 2.5 
30 mm/sec and a step distance of 2.5 mm. The BMC and BMD of the neck, 

Ward's Triangle and trochanteric regions of the proximal femur are 
calculated using a femur software version supplied by Lunar Radiation 
Corporation. The femoral neck region of interest (ROI) is that band about 
1.5 cm wide across the neck of the bone perpendicular to the neutral axis 
35 with the lowest density. Ward's Triangle is defined as a square ROI (about 

1.5 x 1.5 cm) with the lowest density within the proximal femur region. 
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Ward's Triangle is predominantly trabecular bone and contains the least 
amount of bone mineral within the neck region. Carter D.R., et al., 
"Relationship Between Loading History and Femoral Cancellous Bone 
Architecture," J. Biomechanics, Vol. 22, pgs. 231-244, (1989). In the 
5 proximal femur, the ROI is usually the area 1.5 cm wide across the entire 

femoral neck. Additional regions are defined by the software in the lower 
density Ward's Triangle region, and in the region of the greater trochanter. 
Bone loss in the proximal femur begins in the Ward's Triangle region and 
proceeds outward from there. (Brown D., et al., "Mechanical Property 
10 Distributions in the Cancellous Bone of the Human Proximal Femur," Act. 

Orthop. Scand., Vol. 51, pgs. 429-437, (1990). This makes the region an 
early indicator of bone loss, but the higher variance in measuring it, 
compared to the neck region, makes the latter zone, a better discriminator. 
However, Ward's Triangle is the lowest density area at the point where the 
15 neck and greater trochanter meet, a primary hip fracture site. This 

operational ROI may not correspond exactly to the anatomic Ward's 
Triangle region but does provide a repeatable measurement. The width of 
the neck ROI and the size of the Ward's Triangle ROI are actually 
proportional to the measured size of the femoral neck. It has been shown 
20 that the density of the Ward's Triangle area is substantially reduced in hip 

fracture patients compared with age matched controls. (Vose C. P., et al., 
"Femoral Neck Fracturing its Relationship to Radiographic Bone Density," 
J. Gerontol., Vol. 20, pgs. 300-305, (1965). DPA measurements have a 
precision of 1-3% and the scan can be completed in about 20 minutes. Bone 
25 density of vertebrae correlates well with risk for vertebral fracture, and bone 

density of areas in the proximal femur correlates well with risk for hip 
fracture. DPA has a low radiation dose (<10 mrem to skin, and 2 mrera to 
marrow). DPA has approximate sensitivity of spine and hip of 50% and 
53%, respectively, at 95% specificity (i.e. % of fracture cases below 5th 
30 percentile). Moreover, DPA is not subject to systematic errors introduced 

by variable osteoid and variable marrow. 

Bone density can also be measured by radiographic absorptometry. 
Radiographic absorptometry is a method of measuring bone density which 
is well known to those of ordinary skill in the art. Other methods of 
35 measuring bone density include quantitative computed tomography and 

direct measurement of bone density. 
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Direct measurement of bone density can be obtained by measuring 
the bone density in cadavers. Thus, another method of determining the 
correlation between the blood concentrations of certain predetermined blood 
constituents and bone density is to measure the bone density of a set of 
5 cadavers from which blood concentrations of predetermined blood 

constituents are known. The multiple linear regression analysis can then be 
performed and a bone density coefficient relationship can easily be 
determined. 

In correlating the bone density measurements to concentration of 
0 blood constituents, blood concentrations of calcium, phosphate, total 

alkaline phosphatase, an alkaline phosphatase isoenzyme estradiol, and 
progesterone are measured. Liver, bone and intestinal isoenzymes can be 
used. Band 10 alkaline phosphatase isoenzyme. A mathematical 
relationship between the concentrations of blood constituents and bone 
t5 density, as measured by radiographic absorptoraetry or other standard 

method of measuring bone density, is determined by performing a multiple 
regression analysis with the following model: 

Bone density coefficient =* 0 + MCa + b2? + b$E2 + b^H +f>5Alkp + 
20 *>6l-Alkp + £>7(CaCa) + b&(??) + f>9(E2E2) + fclO(Pg + Pg) + 

fcll(AlkpAlkp) + bl2(CaP) + M3(CaE2) + *14(CaPg) +*>i5(CaAlkp) + 
f>16(CaI-Alkp) + fcl7(PE2) + bl8(PPg) + £>19(PAlkp) + i>20(PI-Alkp) + 
&2l(E2Pg) + £>22(E2Alkp) + i>23(E2I-Alkp) + &24(PgAlkp) + fc25(PgI- 
Alkp) + i^6(AlkpI-Alkp) 

25 where: 

Ca = Serum concentration cf calcium 

P = Serum concentration of phosphorus 

E2 = Serum concentration of estradiol 

PG = Serum concentration of progesterone 
30 Alkp = Serum concentration of total alkaline phosphatase 

I-Alkp = Serum concentration of alkaline phosphatase isoenzyme. 

Calculation of the mathematical model utilized the Systat® statistical 

package (Systat: Inc., Evanston, IL). The multiple linear regression 

analysis is an iterative process which calculates the correct coefficients so 
35 that the result of the algorithm using the blood concentrations correlates to a 
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high degree with the result of the bone density measurement by radiographic 
absorptometry. 

For diagnosing osteopenia, the general form of the preferred 
algorithm that is used in the present invention is as follows: 

5 

Bone density coefficient = b 0 + fclCa + fc2P + f>3E2 + btfg +&5Alkp + 
f>6l-Alkp + J>7(CaCa) + &8(PP) + &9(E2E2) + £>lO(Pg + p g) + 
fcn(AlkpAlkp) + £>l2(CaP) + fci3(CaE2) + fcl4(CaPg) +£>i5(CaAlkp) + 
i>16(CaI-Alkp) + fcl7(PE2) + fcl8(PPg) + i>19(PAlkp) + f>20(PI-Alkp) + 
10 fc2l(E2Pg) + £>22(E2Alkp) + £>23(E2I-Alkp) + &24(PgAlkp) + f>25(PgI- 

Alkp) + &26(AlkpI- Alkp) 
where: 

bo is a constant 

Ca = Serum concentration of calcium 
15 P = Serum concentration of phosphorus 

E2 = Serum concentration of estradiol 
PG = Serum concentration of progesterone 
Alkp = Serum concentration of total alkaline phosphatase 
I-Alkp = Serum concentration of alkaline phosphatase isoenzyme. 

20 

It is to be understood that concentration of I- Alkp can be replaced 
with serum concentration of liver, bone or intestinal derived alkaline 
phosphatase isoenzymes. 

In practicing the present invention, the values obtained from the 
25 patient for the six blood constituent concentrations are inserted into the 

algorithm as indicated. The mathematical manipulation is performed. The 
resulting number is called a bone density coefficient and is then placed in the 
severity scale. This results in a bone density probability quotient which 
correlates to a high degree to bone density as measured by radiographic 

30 methods. 

As seen in the algorithm that is considered part of the present 
invention, there are a number of coefficients which are a part of the 
algorithm. It is to be understood that these coefficients can change if the 
bone density coefficient is correlated to a different method of determining 

35 bone density. Thus, if bone density is measured using dual photon 

absorptometry, and one wanted to correlate the bone density coefficient to 



BNSOOCID <WO 9612187A1_L> 



WO 96/12187 



PCT/US95/01379 



49 



10 



the results of the dual photon absorptometry, the overall relationship of the 
tests in the algorithm would be the same or similar as disclosed herein but 
the coefficients could be different 

In its preferred embodiment, the bone density probability quotient is 
assigned to one of the following diagnostic categories: 

Normal to mild osteopenia 
Moderate osteopenia 
Severe osteopenia 

In addition, the menopausal status may be determined based on the 
results of the estradiol and progesterone levels into one of the following 
diagnostic categories: 

15 Probable Prc-menopausal 

Probable Peri-menopausal 
Probable Post-menopausal 

Thus, according to the present invention, using blood concentrations 
20 of certain blood constituents, one can not only diagnose bone density but 

can also obtain an indication of the underlying cause of the osteopenia. 

This invention is further illustrated by the following examples, 
which are not to be construed in any way as imposing limitations upon the 
scope thereof. To the contrary, it is to be clearly understood that resort may 
25 be had to various other embodiments, modifications, and equivalents thereof 

which, after reading the description herein, may suggest themselves to those 
skilled in the art without departing from the spirit of the present invention 
and/or the scope of the appended claims. 

30 Example 8 

In the following example, the bone density coefficient algorithm was 

correlated with bone density as measured by radiographic absorptometry. 

The coefficients were calculated by correlating the blood concentration of the 

indicated blood constituents in 27 patients to the results of bone density in 
35 those patients as measured by radiographic absorptometry. The bone 

density algorithm which results is as follows: 
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0 + (-8.955)(Ca) + (79.370)(P) + (34.076)(E2) + (-9.216)(Pg) + 
(-0.600)(Alkp) + (-57.855X1- Alkp) + (0.926)(CaCa) + (-2.735)(PP) + 
(-0.272XE2E2) + (-0.064)(PgPg) + (0.004)(AlkpAIkp) + (-4.029)(CaP) + 
5 (-3.156)(CaE2) + (1.172)(CaPg) + (-0.031)(CaAlkp) + (8.238)(CaI-Alkp) 

+ (1.906)(PE2) + (-0.138)(PPg) + (-0.018)(PAlkp) + (-6.710)(PI-Alkp) + 
(0.193)(E2Pg) + (-0.048)(E2Alkp) + (-1.825)(E2I-Alkp) 

wherein: 

10 Ca = Serum concentration of calcium 

P = Serum concentration of phosphorus 

E2 = Serum concentration of estradiol 

PG = Serum concentration of progesterone 

Alkp = Serum concentration of total alkaline phosphatase 
15 I-Alkp = Serum concentration of serum lymphycyte-derived alkaline 

phosphatase isoenzyme. 

The bone density coefficient that is obtained from the bone density 
algorithm is assigned to one of three categories: 

20 

Group 1 normal to mild osteopenia 

Group II moderate osteopenia 

Group III severe osteopenia 

Group Ifl <, approximately 78.5 < Group II £ approximately 100 < Group I 
It has been determined that the bone density coefficient calculated according 
to the present invention places patients in the same groups as radiographic 
absorptometry measurements for bone density. Thus, the present invention 
provides a safe, economical and accurate method for diagnosing osteopenia. 

Bone density quotients derived practicing the present invention using 
the six biochemical serum constituents were found to strongly correlate with 
bone density measurements from radiographic absorptometry. It is to be 
understood that the preferred embodiment utilizes all six of the biochemical 
serum constituents listed hereinabove. However, if any five of the serum 
constituents are utilized and then included in the bone density algorithm, a 



25 



30 



35 
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bone density quotient is obtained which does not correlate as well as when 
six tests are used. 

It is contemplated as part of the present invention that any five of the 
biochemical serum constituents can be used to determine the osteopenic state 
5 of a patient For example, if serum calcium is removed from the bone 

density algorithm, then the algorithm predictive value for moderate or severe 
osteopenia is approximately 66% although the predictive value for normal to 
mild osteopenia is 100%. If progesterone is removed from the bone density 
algorithm, then the algorithm predictive value for moderate or severe 

10 osteopenia of approximately 75% and a predictive value for normal to mild 

osteopenia of approximately 94%. If alkaline phosphatase is removed from 
the bone density algorithm, then the algorithm predictive value for moderate 
or severe osteopenia of approximately 83% and a predictive value for 
normal to mild osteopenia of approximately 85%. If estradiol is removed 

15 from the bone density algorithm, then the algorithm predictive value for 

moderate or severe osteopenia of approximately 67% and a predictive value 
for normal to mild osteopenia of approximately 85%. However, it has been 
determined that when all six tests are included in the bone density algorithm, 
the predictive value for moderate and severe osteopenia is 100% and the 

20 predictive value for normal to mild osteopenia is 100%. All of the predictive 

values may vary slightly with different patient populations. 

Example 9 

The following serum constituents are measured in a 37 year old 
25 woman: 



Calcium 9.5mg/dL 

Phosphate 5.0mg/dL 

Estradiol ; 2pg/mL 

30 Progesterone 0.2ng/mL 

Alkaline Phosphatase 80 U/L 

Alkaline phosphatase isoenzyme 0 (negative) 



The above blood chemistry values were used to calculate a bone 
density coefficient using the following algorithm: 
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0 + (-8.955)(Ca) + (79.370)(P) + (34.076)(E2) + (-9.216)(Pg) + 
(-0.600)(Alkp) + (-57.855)(I-Alkp) + (0.926)(CaCa) + (-2.735)(PP) + 
(-0.272)(E2E2) + (-0.064)(PgPg) + (0.004)(AlkpAlkp) + (-4.029)(CaP) 4- 
(-3.156)(CaE2) + (L172)(CaPg) + (-0.031)(CaAlkp) + (8.238)(CaI-Alkp) 
+ (1.906)(PE2) + (-0.138)(PPg) + (-0.018)(PAlkp) + (-6.710)(PI-Alkp) + 
(0.193)(E2Pg) + (-0.048)(E2Alkp) + (-1.825)(E2I-Alkp) 



where: 

Ca = Serum concentration of calcium 

P = Serum concentration of phosphorus 

E2 = Serum concentration of estradiol 

PG = Serum concentration of progesterone 

Alkp = Serum concentration of total alkaline phosphatase 

I-Alkp = Serum concentration of serum lymphycyte-derived alkaline 

phosphatase isoenzyme. 



A bone density coefficient of 101.2 is calculated for this patient. 
This coefficient falls within Group I in the following severity scale: 

Group in < 78.5 < Group H < 100 < Group I 



This patient has normal to mild osteopenia. When the patient's bone 
density is measured by radiographic absorptoraetry, she is found to have 
normal bone density. 

Example 1ft 

The following serum constituents are measured in a 47 year old 
woman: 



15 



20 



Group I 
Group II 
Group HI 



normal to mild osteopenia 
moderate osteopenia 
severe osteopenia 
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Calcium ■ 9.5mg/dL 

Phosphate 5.0 mg/dL 

Estradiol 2pg/raL 

Progesterone 0.2 ng/ml 

5 Alkaline Phosphatase 80 U/L 

Alkaline phosphatase isoenzyme 1 (positive) 



The values are inserted into the bone density algorithm recited in 
Example I and a bone density coefficient of 84.4 is calculated. This falls in 
Group II, indicating a moderate osteopenia. When the patient's bone 
density was measured by radiographic absorptometry, she was found to be 
moderately osteopenia 

Example 11 

The following serum constituents are measured in a 47 year old 
woman who has been diagnosed with breast cancer: 



Calcium 9.0rag/dL 

Phosphate 4.0 mg/dL 

20 Estradiol Opg/mL 

Progesterone.. 0 ng/mL 

Alkaline Phosphatase 60 U/L 

Alkaline phosphatase isoenzyme 1 (positive) 



23 The values are inserted into the bone density algorithm recited in 

Example I and a bone density coefficient of 69.9 is calculated. This falls in 
Group m indicating a severe osteopenia insured by another method. 

Example 12 

30 The following example shows the correlation between the method of 

the present invention and Dual Photon Absorptometry (DP A) as a tool for 
diagnosing and determining the severity of osteoporosis in a particular 
individual. 200 female subjects representing a cross section of ages and 
menopausal status were evaluated. These specific age groups were 

35 represented as follows: 
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Ag£ Number of PflfjflUfi 

25-35 26 

36-45 25 

46-55 33 

5 56-65 30 

66-75 34 

76-85 25 

85 and over 26 



10 Each subject completed an Osteoporosis Data Questionnaire Form, 

which included detailed information on personal history, family history, 
gynecologic history, medical history, surgical history, and drug history, 
with special references to hormone treatment. 

Three tubes of blood were collected from each subject via 

15 venipuncture. One tube was sent to Barnhill-MetPath Laboratories in 

Savannah, Georgia for analysis of serum calcium, phosphorous, and total 
alkaline phosphatase. The second tube was sent to Horus Therapeutics, 
Inc. in Savannah, Georgia for determination of lymphocyte alkaline 
phosphatase utilizing the Isoelectric Focusing Method. The third tube was 

20 sent to MetPath Laboratories in Teterboro, New Jersey for evaluation of 

estradiol and progesterone by immunoassay. Data from all three sources 
were then subjected to statistical analysis. At the time of the analysis of 
data, each subjects chart was critically reviewed for relevant clinical 
information. 

25 The Bone Mineral Density (BMD) of each subject was also 

measured by Dual-Photon Absorptometry using 15 3 Gadolinium as a source 
(using Lunar DP3, Lunar Radiation Corporation, Madison, Wisconsin). 
DP A measurements were performed on the Lumbar Spine (L1-L4) and the 
hip, specifically the femoral neck, Ward's Triangle, and the trochanteric 

30 regions. BMD was measured in gra/cm 2 for each region noted, and the 

fracture risk was determined. 

The BMD measurements were then age matched, and adjusted for 
sex, age, ethnic group, and weight. From the individual values obtained, 
average BMD for hip and spine was calculated. The overall average BMD 

35 for both hip and spine was also computed. 
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A mathematical algorithm correlating the results using the present 
invention with DPA was determined with the resultant data is shown in 
Figure 11. 

Figure 12 illustrates the discrepancies between the DPA 
measurements and the measurements obtained using the method of the 
present invention for individual subjects. The amount of the discrepancy 
(the residual) for any patient is the vertical distance between the points for 
that patient on the DPA and the present invention lines. 

Figure 13 illustrates the appropriateness of using the present 
invention as a prediction of DPA measurements. Note the absence of any 
tendency for points to be consistently above or consistently below the 
diagonal straight line, as greater DPA measurements are considered. 
Additionally, there is an absence of any tendency for points along any 
vertical line to be more widespread than are points along any other vertical 
line. Finally, note that the generally elongated shape formed by the points is 
consistent with a strong linear correlation between DPA and the present 
invention. 

Figure 14 demonstrates the normality of values predicted by the 
present invention (the horizontal scale). The vertical scale shows the of the 
normalized scale (z) scores values determined according to the present 
invention as would be expected from the assumption of the present 
invention values forming a normally distributed collection of numbers. 

The closer the dots are to the straight line, the greater is the 
consistency between what is observed about the present invention values 
and what is expected about them from the assumption of their forming a 
normally distributed collection of numbers. 

Figure 15 demonstrates the normality of residuals (the difference 
between the observed DPA measurement and the value predicated by present 
invention), as evidenced by the arrangements of points along the diagonal 
straight line. The horizontal scale shows values of residuals The vertical 
scale shows the residuals' normalized scale (z) scores as would be expected 
from the assumption that the residuals form a normally distributed collection 
of numbers. 

The closer the dots are to the straight line, the greater is the 
consistency between what is observed about the residuals and what is 
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expected about them from the assumption of their forming a normally 
distributed collection of numbers. 

Figure 16 demonstrates that the variance of residuals stays the same 
regardless of the value calculated according to the present invention involved 
5 in the residuals' computation. This homogeneity of variance appears in the 

uniform density of the dots in a rectangle formed by residuals between -0.17 
and +0.17 and by the present invention between 40.4 and +1.0. Moreover, 
there is no discemable tendency for points along any vertical line to be more 
widespread than arc points along any other vertical line. 
10 Figure 17 demonstrates the independence of residuals and the values 

predicted by the present invention. The horizontal scales shows the values 
predicted by the present invention. The vertical scale shows the residual's 
studentized scale (t) scores. 

When there is no relationship between values predicted by the 
present invention and the residuals involving those predicted values (that is, 
when there is independence between residuals and predicted present 
invention values) the following criteria should be met: 

1 . There should be no discemable pattern (such as a straight or 
curved line) to the dots. 

2. About half the dots should appear in, the top half of the 
graph. 

Figure 18 demonstrates that virtually all observed DPA 
measurements are satisfactorily predicted by the method of the present 
invention. The horizontal scale shows values predicted by the method of the 
present invention. The vertical scale shows Cook's distances between 
observed DPA measurements and the predictions of those measurements by 
the method of the present invention Thus, the closer the dots are to a height 
of Cook = 0, the more accurate is the method of the present invention. 
Because Cook's distances tend to form a collection of numbers fitting an F- 
statistical distribution, the two outline points at heights between 10 and 20 
are not thought to be problematic in view of the large size of the sample. 



20 



25 



BNSDOCID: <WO 9612187A1J > 



WO 96/12187 



PCTAJS95/01J79 



57 

Example 1? 

The following example describes the training of a neural network to 
prognose prostate cancer. 

As shown in Fig. 19, a total of 52 samples were divided into 2 

5 groups, a training set and a generalization testing set. The training set 

contained 40 samples (28 stable and 12 progressing) and the generalization 
testing set contained 12 samples (9 stable and 3 progressing). 

The initial network architecture was selected based on the level of 
complexity of the classification task. A multi-layer feedforward network 

10 was used. Selection of the initial architecture involved the selection of the 

number of hidden layers and the number of neurons in each hidden layer. 
Several trial iterations were performed to determine and adequate 
configuration that showed good results on both the training sample set and 
the generalization test sample set. The present network had one hidden 

15 layer, having nine neurons, and two output neurons. 

Initially, connection weights among the neurons were randomly set 
The neural network had five input neurons, corresponding to five input 
variables significant for prostate cancer: TPS, PSA, PAP, CEA, and 
testosterone. The training data is shown in Fig. 19. During training, the 

20 five input variables for each patient were first linearly scaled into the 

continuos range between 0.0 and 1.0. The resultant five numbers were then 
presented as an input vector to the input neurons of the artificial neural 
network. 

For each of the input vectors, the network generated an output based 
25 on the connection weights among the network neurons. The output can be a 

single value or a vector of numbers, depending on the number of output 
neurons used. The network used had two output neurons. The outputs of 
the two neurons were processed by the following mathematical equation to 
yield a single diagnostic index: 

30 

Index = f ANN2-ANNH + 0.5 
2 

Each neuron in the network participated in the output calculation by passing 
35 the sum of all inputs to the neuron through a non-linear s-shaped function 

(often a ligistic function) and sending the result to each and everyone of the 
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neurons in the following adjacent layer. The generated output or each 
output neuron was compared to the desired "target" output. A value of 0.1 
corresponded to a diagnosis of stable and an output of 0.9 corresponded to a 
diagnosis of progressing. The difference was used to calculate an error terra 
to guide the training algorithm, i.e., the back propagation algorithm, in the 
adjustment of network connection weights in an attempt to reduce the 
differences between network outputs and target values over the training 
sample set. 

After training, the neural network correctly classified 100% of the 
samples. 

When presented with the generalization test results, the trained 
neural network correctly identified 100% of the stable samples and 66% of 
the samples where the disease was progressing. 

It should be understood, of course, that the foregoing relates only to 
preferred embodiments of the present invention and that numerous 
modifications or alterations may be made therein without departing from the 
spirit and the scope of the invention as set forth in the appended claims. 
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Claims 

What is claimed is: 

1 . A method for diagnosing a disease in a human or animal 
comprising the steps of: 

measuring the concentrations of a predetermined set of biomarkers 
known to be associated with the disease from a biological fluid from 
the human or animal; 

scaling the digitized values of the biomarkers concentrations; and 

introducing the scaled values to a trained neural network means, 
whereby the output values from the neural network means tend 
toward the upper value when the human or animal has the disease 
and the output values tend toward the lower value when the human 
or animal does not have the disease. 

2 . An apparatus for diagnosing a disease in a human or animal 
comprising: 

means for digitizing the concentrations of a predetermined set of 
biomarkers known to be associated with the disease from a 
biological fluid from the human or animal; 

a means for scaling the digitized values; and ' 

a trained neural network means coupled to the digitizing and scaling 
means for generating network output values between an upper and 
lower value, whereby the output values tend toward the upper value 
when the human or animal has the disease and the output values tend 
toward the lower value when the human or animal does not have the 
disease. 
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3. A method for diagnosing a disease in a human or animal 
comprising the steps of: 

5 measuring the concentrations of a predetermined set of biomarkers 

known to be associated with the disease from a biological fluid from 
the human or animal; 

scaling the digitized values of the biomarkers concentrations; 

10 

introducing the scaled values to a first trained neural network means, 
whereby the output values from the first neural network means tend 
toward the upper value when the human or animal has the disease 
and the output values tend toward the lower value when the human 
15 or animal does not have the disease; and 

introducing the output value from the first neural network and a 
second set of predetermined biomarkers to a second trained neural 
network means, whereby the output values from the second neural 
20 network means tend toward the upper value when the human or 

animal has the disease and the output values from the second neural 
network means tend toward the lower value when the human or 
animal does not have the disease. 



25 



4. The method of Claim 3 wherein the second set of 
predetermined biomarkers includes patient biographical information. 



5. The method of Claim 3 wherein the second set of 
predetermined biomarkers are the same as the first set of predetermined 
30 biomarkers. 
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Neural Network Diagnostic System 
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105 


60 


9.1 


3.1 


29 


0.3 


88 


10.7 


46.023 


9.631 


106 


77 


9.1 


4 


15 


0.4 


66 


5.8 


68.636 


0.603 


107 


46 


8.8 


4 


14 


0.2 


44 


o 


72.500 


0.784 


108 


52 


8.5 


3.2 


14 


0.2 


88 


6.1 


43.409 


0.663 


109 


52 


8.6 


3.7 


10 


0.1 


62 


3.8 


41.452 


0.700 


110 


79 


9.1 


3.8 


13 


0.2 


78 


o 


71.795 


0.587 


111 


66 


9.5 


4.4 


9 


0.1 


74 


o 


64.189 


0.652 


112 


40 


9.2 


3,1 


3? 


9.7 


51 


6.8 


41,961 


0-758 


113 


82 


9.7 


3.7 


5 


9.1 


70 


5,8 


72,714 


o? 70 


114 


66 


10 


4.5 


9 


9.1 


66 


0 


67,879 


0.662 H 


115 


62 


9.3 


3 


34 


9,1 


63 


0 


71.587 


0.6?0 


116 


60 


8.7 


3.9 


14 


9.1 


56 


11,7 


59-464 


0',6frQ 


117 


81 


8.8 


3.6 


17 


0.2 


96 


7,3 


M.792 


9,>tt 


ft 118 


75 


9.2 


3.3 


11 


0.2 


72 


9 


61.528 


0.515 


ft 119 


66 


8.9 


4.5 


9 


0.2 


81 


o 


45.062 


9,65 1 


H 120 


78 


8.4 


3.9 


7 


0.1 


71 


8.5 


62.676 


0.501 


121 


74 


8.8 


2.8 


0 


0.1 


62 


o 


62.742 


0.502 


122 


81 


8.8 


3.2 


10 


0.1 


61 


6.8 


28.933 


0.581 


123 


74 


10.4 


3.4 


11 


0.2 


69 


0 


48.551 


0.600 


124 


82 


8.8 


3.4 


17 


0.2 


79 


o 


61.772 


0.560 


125 


44 


8.6 


3.2 


208 


0.3 


120 


o 


70.417 


0.728 


126 


69 


8.5 


4.5 


18 


0.1 


78 


5.9 


57.179 


0.634 


127 


75 


8.7 


3.8 


9 


0.2 


69 


0 


66.957 


0.605 


12? 


74 


9 


3.8 


8 


9,1 


89 


4 


72.375 


0.599 


139. 


89 


3,7 


2.9 


32 


9,2 


139 


9. 


68.912- 


ft S")^ 1 


131 


77 


9.5 


4.2 


41 


0.1 


74 


6.8 


69.054 


960M 


132 


72 


8.9 


3,5 


24 


9,3 


58 


o 


43.621 


0.620 1 


123 


6? 


9.2 


3,1 


9 


0.2 


94 


8,4 


49.468 


9,59 7 1 


liL 


71 


9 


6.3 


5 


9.1 


79 


4,6 


$1,571 


9,646 


135 


ft 


9.2 


3.6 


11 


9.1 


84 


o 


72.857 


9,*!6_ 


1?6 


85 


9 


3,9 


7 


9.1 


71 


12,6 


M.366 


9-571 


137 


43 


9,2 


3.8 


149 


0.2 


46 


9 


68,261 


0.777 


138 


49 


9.2 


2.9 


73 


3.3 


42 


0 


57.619 


0.696 


139 


41 


8.6 


4.1 


117 


15.7 


93 


o 


63,656 


0.716 


J 140_ 


67 


9 


3.4 


93 


0.1 


59 


o 


73.051 


9,646 


\ 141 


59 


9.8 


2.1 


132 


0.2 


45 


3.5 


42.000 


0.702 


B tti 


70 




2^7 


XL 


0.1 


f£L 


\ 0 


72.985 


0.600 



SUBSTITUTE SHEET (RUl£ 26) 

BNSOOCID: <WO 9612187A1 J. > 



WO 96/12187 



8/24 



PCT/US95/01379 



■ 

PATIENT 


AGE 


CAL. 


PHOS 


1 1 wo 


' 11 
PROG_ 


TAI D 




pi rvpR 


SCORE 


143 


66 


9.4 


3.3 


1 A 




/ 1 


ft 


71 690 


0 625 


144 


52 


9 


3.5 








R 


35 156 


0.706 


* A C 

145 


82 


8*2- 




1 c 


n i 




o 


49 063 


0.572 


146 


SLA 

o4 






140 


ft S 


47 


8 7 


58.936 


0.675 


« A 1 

147 


65 




.., 4i6- 


o 


ft ** 


78 


o 


41 154 


0.661 


% AO 

148 


?J_ 


7.4.. 




1 A 


ft ? 


96 


5 7 


69.583 


0.683 


t A f\ 

149 


A Q 


2i2_ 


•♦.1 




ft ^ 


143 


8 7 


34 895 


0.708 


150 


75 


10 


3.9 


19 


o,2 


8? 


o 


39,545 




151 


57 


9,4 


3,9 


10 


0,2 


54 


3.4 


51.111 




151. 


47 


9,4. 


4 


42 


0,3 


65 


8,6 


5I.M6 


0.739 


n .i53. 


70 


8.8 


3.5 


6 


0.1 


65 


5 


37.231 


0.621 


154 


59 


?.2 


4 


9 


0.2 


49 


4.4 


47.143 


0.690 


155 


60 


9.2 


4.1 


12 


0.1 


56 


2.8 


60.357 


0.680 


156 


49 


?,1 


3.6 


59 


0,2 


35 


0 


54.571 


0.762 


157 


59 


9,5 


4.1 


7 


0.1 


92 


6.5 


45,975 


0.660 


158 


66 


9,2 


3.8 




0.2 


46 


0 


40.000 


0.660 


159 


57 


9.2 


4.3 


16 


0.2 


81 


5,4 


53.704 


0.679 


160 


41 


10 


4.6 


306 


1.2 


46 


o 


75.217 


0.763 


161 


3? 


9.5 


3,6 


4 


0,2 


6? 


5.1 


69,26? 


0.758 


162 


35 


9,5 


3.4 


36 


0.4 


57 


3.3 


64.561 


0.785 




SUBSTITUTE SHEET (RULE 26) 

BNSOOCID: <WO 9612187A1_I_> 



WO 96/12187 



9/24 



PCTAJS95/01379 



Testing Data for Osteoporosis 
Diagnostic System 



II PATIPNT 1 AGE 


CAL 1 


PHQ5 


ETWQ 1 


PROG 1 TALP 1 


T1NTES 


PLIVER 


SCORE 


m 


52 


9,4 


4.2 


-1 


0.3 1 


61 


o 


46 


P.722 


II 128 1 


27 


9,8 - 


4.4 1 


-1 


0.4 1 


58 


o 


43 


0.862 


H 121 1 


39- 


9,1 | 


3,9 1 


15$ . 


11.5 1 


44 


0 


59 


0.820 


ll m l 


26 - 


9,5 1 


3,8 1 , 


41 


0.2 1 


98 


0 


$9 


0,790 




201 


30 


9,1 1 


3.5 1 


123 


3.3 1 


90 


o 


52 


0.814 




206 1 


39. 


9.2 1 


2.6 


50 . 


1 S 


74 


89 


42 


0.729 




213 1 


46 


9.7 1 


3.4 1 


37 


0 7 1 


90 


o 


36 


0.708 




220 1 


4$ 


8.8 1 


3.4 | 


118 


4 3 1 


90 


9 


63 


0.996 




•m 1 


49 


9.7 I 


3.2 1 


72 


0 4 1 


78 


8.6 


38 


0.717 




7*4 38 


9.5 1 


3.2 1 


186 


12.6 1 


50 


o 


46 


0.803 




243 


22 


9.3 1 


4 .4 I 


« 


0.3 1 


58 


0 


40 


0.857 




244 1 


53 


9.6 1 


3,7 




0.3 1 


58 


0 


49 


0.725 




245 


39 


9.5 1 


3.2 1 


-J 


0.1 1 


64 


o 


53 


0.800 




747 1 


28 


9,1 | 


3,$ 


-1 


0.3 1 


58 


0 


50 


0.821 


1 249 22 


9,2 | 


3,5 | 


42 


0.2 1 


9? 


11.1 


38 


0.782 


753 1 24 


9.2 1. 


3,4 1 


-1 


0,3 1 


« 


o 


63 


0.797 




254 30 


9.el 


4 1 


142.. 


10 1 1 


45 


o 


21 


0,894 




256 1 Li. 


9.1 


3,8 1 


219 


0.2 1 


5? 


6,1 


32 


0.871 




258 1 30 


9,3 


3,4 1 




0.$ 1 


68 


o 


33 


0.829 




259 1 21 


9,2.1 


4.4 1 


44 


0 1 1 


108 


o 


54 


0.805 




760 1 23 


9.1 


3.8 


-1 


0 % 1 


62 


o 


61 


0818 




7A1 


73 


9.4 1 


3.4 J 


-I 


0 2 1 


85 


o 


32 


0.795 




265 24 


9.4 1 


3.2 1 


191 


0 2 I 


71 


0 


30 


0.829 




266 1 26 


9 1 


3.9 1 


91 


0 3 1 


51 


5.5 


47 


0.869 




268 1 20 


9.3 1 


4.3 1 


62 


0.2 i 


61 


0 


60 


0.857 




269 1 L2_ 


8.9 1 


3.9 | 


53 


0.8 1 


68 


o 


42 


0.837 




271 I 20 


8.9 1 


4.1 1 


277 


18.9 1 


65 


o 


40 


0.869 




777 1 74 


9 1 


3.3 1 


58 


0.4 1 


66 


o 


49 


n 01C 

0.815 




774 


19 


9.7 1 


3,8 1 


44 


0 3 1 


120 


o 


70 


0.776 




276 25 


9.2 1 


3.9 1 


-1 


0.9 


65 


5 


27 


0.837 




770 1 71 


10 i 


4,$ 


115 


0.5 1 


106 


o 


70 


0.834 




280 1 
281 

284 


21. 
20 


9.5 . 
9.1 . 
9 


A 1 

4.1_ 
3,6 

I 3,3 


7ft7 
1 'I 

37$ 


0.9 1 

0.3 

0.6 1 


$8 

47 

42 


o 

o 
o 


57 
53 
37 


0.856 

0.82J 
0.855 




285 

286 

7R7 


| 22. 

1 31 

27 

1 77 


9.6. 
9 

9_ 


1 3 - 7 
1 3,2 

1 31 


208 
350 

45 


1.2 1 
1.5 I 

0,2 | 


to 

65 
66 


o 


62 
60 
39 


0.842 

. 0,?19. 
0.805 




290 23 


9.3. 


4,2 


1 57 


0,3 1 


104 




47 


0.817 




707 1 74 


9.4 


1 3.7 


1 35 


0.3 1 


ft 




43 


0.829 




296 


2| 


. 9.2 


1 32 




0.3 


111 


o 


55 


0.767 




7QR 


1 76 


9.4 


1 3.2 


248 


0.8 


57 




35 


0.842 




321 39 


8.9 


1 3.6 


J 84 


0.2 


ft 


o 


42 


0.785 




326 1 22 


9.4 


1 3.6 


4J 


0.4 


50 


o 


37 


0.847 




331 1 48 


. 9.7 


1 31 


1 47 


0.3 


84 




58 


0677 




332 1 51 


. 9.8 


1 3.5 


- 1 


0.3 


64 




40 


0707 




XVK 1 5* 


9.2 


1 3.3 


1 166 


1.2 


63 




56 


0.705 




340 


1 42 


n, 


1 2.8 


L™5J 


9,?. 


L 


0 


64 


9 728 



<WO 9612187A1 I > 
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SUBSTITUTE SHEET (RULE 26) 



WO 96/12187 
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Arc 

Aire 






PTWO 




TA1 P 




PI TVFR 






2y_. 


0 ft 

y .o 


1 R 


7ft 


0 1 


ft7 
oz 


n 

v 


71 

/ 1 


0 710 

v- / jU 


JHQ 


44 


0 7 


1 4 


1 1 1 
III 


A Q 


77 

l z 


17 A 


SO 


O Q77 


\A1 


41 


Q S 


3 7 


.1 


0 7 


77 


o 


ftQ 


0 71S 


3S0 

J J 7 


54 


0 3 




OS 


0 2 


77 


o 


6Q 


0 603 


IftO 


34 


0 4 


3 ft 


67 


0 4 


1(\ 


n 


63 


0 701 


1ft4 


40 


10 


4 7 


104 


0 1 


ft7 
oz 


Q 


S7 


0 fti 1 


IftS 


4ft 


o 

7 


1 S 


07 

y / 


ft 7 


<4 


A 

1L 


ft4 
VC- 


O ftOft 


17s 

j / j 


si 

22- 


0 1 


1 0 

1x7 


7* 

L2- 


O 7 


ft7 


fi 

U- 


ftft 

SSL. 


O 71ft 


1fl7 
38Z 


fLa_ 


ft 7 
o. / 


A ft 


3ft 


O 7 

ILL. 


fti 


A 

LL 


ftft 
oo 


O 70ft 


Iftft 


74 
Zh 


O 1 


1 fi 
3.0 


.i 

*4 


O 1 

V.J 


ftft 


A 
U 


ftS 


O ftl7 


3oy 


77 
6*- 


O 1 

y_j 


4 1 

H.J 


" 1 


O 1 


ftft 

DO 


A 


ftS 

— _ — Oj 


O ftll 


101 

3yi_ 


17 


O 7 


1 1 
2l2- 


.1 


0 4 


S7 


A 

V 


ft4 


0 707 

— \). tyi 


101 


7Q 

_ly 


O ft 

y.o 


1 S 


IO 

31/ 


0 s 


47 
H / 


A 

v 


fti 


O ft7ft 


^73- 


77 


ft Q 

o.y 


4 1 

H.I 


.1 
*1 


O 4 
U.H 


4S 


A 

v 


Sft 

20 


O ft4ft 


4A7 
4UZ 


39, 


ft 7 


1 S 

J.J 


11 

j j 


O ft 


41 


A 


4ft 


O 7ftft 


401 
hU3 


17 
3Z 


O 7 

y.z 


1 4 

3.H 


i 

-I 


0 1 

v 1 J 


71 
Z 1 


A 

V 


74 

£=- 


O ft47 

U.oHZ 


Add 
huh 


40 
22- 


0 7 

y .z 


4 S 

H.J 




0 1 


100 


A 

li- 


SI 

22- 


O 71ft 


407 

H\J / 


S7 


0 4 


1 1 

_ 3,3 


1S7 

LiZ 


0 4 

v.H 


71 

Li- 


ft 


S7 
2Z 


0 704 

\J . / vH 


408 


7d 


Q 


3 7 




0 ft 


ft? 


o 


40 


0 810 

U.017 


^NV 


Aft 


0 7 


4 1 

H , 1 


ftO 


7 1 


S7 


o 


Sft 


0 608 


417 


S7 


0 7 


1 ft 


1 


0 7 


ftS 


o 


70 


0 704 


41ft 


44 

hh 


0 7 


1 4 


147 

IHi 


0 7 


80 


o 


70 


0 602 


4 i'- 


40 


O S 

y.j 


1 4 


> 


0 7 


7S 
li- 


U- 


47 

H / 


0 70ft 


4 15 


AQ 


IO 


A 7 

H.Z 


401 


1 1 




A 

y- 


7S 


O ft7ft 




oft 


O 1 




t 


0 A 


07 
0Z 


A 

Li. 


ftS 

Oj 


O 7ft7 


417 
HZZ 


ift 


O 7 
7i*_ 


1 1 


lOft 

iyo 


O 7 

u.z 


7ft 
/O 


A 

u 


74 

^H 


O 7RO 


SIO 

j3y 


4ft 

ho 


O 1 

y*3_ 


4 * 

H.J 


701 

ly 3 


O 7 

v./ 


4ft 


A 

U 


11 


0 813 

V,0 1 3 




7ft 


0 7 


1 
3 


Oft 


0 1 


ftS 




o 

V 


47 


0 82S 


S40 

D^y 


7ft 


O S 

y,j 


1 O 

3.y 


SO 


0 4 

vF.H 


4ft 

HO 


7 7 


40 


0 862 


SftO 


ii 


fi 0 


1 4 


i ift 


0 1 
lU— 


70 

iy 


7 4 


61 


0 701 


Sft7 


71 
—jfc2 . 


O 4 

-Li- 


1 7 
3.Z 


i 


0 4 


Sft 

2&_ 


n 

V 


JO 


0 81S 


579 


7A 


ft 


1 4 

3.4 


48 


n < 

U.5 


ftO 


n 





O 707 


<QQ 

588 


36 


A 1 

9.3 


1 * 

3. j 


i 


O 1 

ILL 


4ft 
48 


n 
v 


«i 

31 


O 701 




7Q 

Zo 


O 7 

y. / 


1 7 

3. / 


1ft 
3D 


o < 

U.J 


S7 

jZ 


A 
V 




0 8^8 


593 


35 


9 


3.7 


108 


0.3 


69 


9.6 


41 


0.811 


597 


3? 


9.3 


3,5 


-1 


0.3 


53 


0 


71 


0,801 


m 


49 


9,3 


3,6 


180 


4,3 


$9 


0 


46 


0.693 


2520 


29 


9.2 


4.4 


52 


Q,2 


49 


0 


41 


0,«77 


3000 


28 


9.4 


3,7 


4? 


0.4 




0 


53 


0.828 


3020 


25 


9,$ 


4 


164 


16,8 




o 


41 


0,880 


100 


52 


5.7 


3,1 


....37J. 


0,3 


59 


0 


54 


0.703 


113 


53 


9.8 


3.5 


3 


0.2 


63 


o 


67 


0.690 


126 


42 


9 


4 


99 


16,4 


83 


14.4 


36 


0 984 


130 


50 


9.3 


3.5 


336 


0.6 


40 


0 


25 


0 754 


139 


74 


8.9 


3.7 


209 


0.4 


35 


5.9 


36 


0W8 


160 


W 


9.5 


3.9 


81 


0.3 


47 


5.2 


53 


0'04 


176 


-_56 


? 


2.9 


115 


0.2 


64 


o 


58 


0 h77 




SUBSTITUTE SHEET (RULE 26) 
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WO 96/12187 



11/24 
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., ■ _ ... 
PATIENT 


AGE 


CAL_._ 


PHQS 


ETWO 


PROG 


TALP 


..TINTES 


PLIVER 


SCORE 


177 1 


?9 I 


9,4 


3.1 


71 


0.2 


82 


o 


45 


A T< 1 

0.751 


208 1 


42 1 


9.1 . 


2.4 


177 


15.6 


54 


o 


•?A 

70 


,,-9.757 


211 1 


50 1 


9.5 


3.1 


69 


0.2 


77 


O 


JCA 


U.oVO 


214 1 


25 1 


9,2 . 


3.7 


-1 


0.9 


43 


0 


*>o 
ly 


0.852 


25J_L 


24 


9,5 


3.9 


-1 


0.1 


61 


0 


25 


A QAA 

0.844 


273 


30 


9.4 


3 


99 


0.2 


156 


8.3 


V* 


0.776 


299 1 


27 1 


9.2 


3.2 


41 


0.2 


77 


0 


58 


A "TOO 

0.788 


309 1 


35. . 


9 


3.9 


55 


0.3 


90 


8.7 


62 


9.7Z2- 


316 1 


41 


9.4 


3.5 


75 


7.3 


A A 

44 


0 


I / 


0.733 


320 1 


52 1 


9.3 


3.4 


57 


0.2 


58 


r\ 
O 


<< 

jj 


A TAT 

o. /u/ 


355 | 


54 


9.1 


3-6 


93 


0.4 


55 


0 


40 


A 1 1 0 

0.718 


358 1 


54 1 


9,5 


3.6 


38 


0.3 


71 


0 


54 


A iCOO 

0.688 


361 1 


4? 


9 


3.8 


31 


0.3 


74 


9.5 


50 


s\ i A A 

0.744 


371 1 


4? 


9.5 


3.8 


-1 


0.2 


72 


10.3 


42 


0.698 


377 1 


52 1 


9.$ 


3-6 


129 


0.2 


68 


0 


53 


0.717 


385 1 


55 


9 


3.3 


-1 


0.2 


97 


o 


73 


A dlA£. 

0.646 


398 I 


5? 


9,3 


4.6 


30 


0.2 


97 


1Q.? 


28_ 


A JlA^T 


4Q6 1 


2$ 


9.6 


3.4 


35 


0.3 


$2 


0 


"TA 

70 


A QAA 

O.oOo 


414 1 


5? 


9.3 


3-9 


56 


0.4 




o 


72 


A "7 I A 

0.710 


421 


31 


9 


3,5 


72 
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• * Program 


: Calc. prg 


Author 


: Terminal solutions (JJK) 


•* Date 


: 14 August 1992 


** Purpose 


: Quios calculating program 



set decimal to 16 
set fixed on 

•**•** Set Coefficient information to valiabtes 
store 0.000000000 to outnumb 
store - 1009.33 to constant 
store - 0.014786248 to c _age 

store 424.462000000 to c_ cal 
store - 66.3468 to c_cal2 

store - 0.025761569 
store - 0.019715699 
store 0.027232306 
store - 0.000000263 
store - 12.0319 
store 7.499334732 
store - 2.220366330 
store 0.313008013 
store - 0.016909117 
store 0.001369629 
store - 0.000428831 
store 0.001798351 
store 0.001841052 
store 0.007649946 
store 0.002153604 
store 4.606320653 
store - 0.119876531 
store - 0.011786744 
store - 0.000170224 
store 0.000000706 
outnumb= (constant) 
outnumb=outnumb + (c_age * age) 
outnumb=outnumb + (c_ cal * cal) 
outtiumb=outnumb + (c_ cal2 * (cal * cal)) 
outnumb=outmimb + (c_ csage • (cos(age))) 
outnumb=outnumb + (c_ csphos2 * (cos(phos * phos))) 
outnumb=outnumb + (c_ cstintcs * (cos(tintes))) 
outnumb^outnumb + (c_ ctwoimx2 * (etwoimx • etwoirax)) 
outnumb=outnumb + (c_ phos * phos) 
outnumb=outnumb + (c_ phos2 * (phos * phos)) 
outnumb=outnumb + (c_ phos3 * (phos * phos. * phos)) 
outnumb^outnumb + (c_ phos4 * (phos * phos * phos * phos)) 
outnumb=oumumb + (c_ phos5 * (phos • phos * phos * phos * phos)) 
putnumb^outnumb + (c_ prog2 * (prog * prog)) 
oumumb^outnumb + (c_ progpliv * (prog * pliv)) 
outnumb^outnumb + (c_ progtint * (prog • tintes)) 
outnumb=outnumb + (c_ metwoimx • (tan (etwoimx))) 
outnumb=outnumb + (c_ tntintes * (tan(tintes))) 
outnumb=outnumb + (c_ agephos * (age * phos)) 
outnumb=outnumb + (c_ cal3 • (cal • cal * cal )) 
outnumb=outnumb + (c_ cal4 • (cal * cal • cal * cal)) 
outnumb=outnumb + (c_ talp * talp) 
outnumb=outnumb + (c_ talp2 * (talp * talp)) 
outnurob=outmimb + (c_ talp3 * (talp * talp * talp)) 
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to c_csage 
to c_csphos2 
to c_cstintes 
to c_etwoimx2 
to c_ phos 
to c_ phos2 
to c_ phos3 
to c_ phos4 
to c_ phos5 
to c_ prog2 
to c_ progpliv 



to c_ progtint M/^MTri "/ ~] 

# / 



to c_ tnetwoimx 
to c_ tntintes 
to c_ agephos 
to c_ cal3 
to c_ cal4 
to c_ talp 
to c_talp2 
to c_talp3 
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N 4 

Nl 8 

N2 1 

f(x) L 



i = 



J = 



O N-l 4 input nodes 
0 Nl I 8 neurons in first hidden layer 
1 output neurons 
Logistic neuron output function 



1+e 





100 




0 


I max. = 


25 


Imin. = 


10 




40 




0 




40 




0 



I range. = I max. -Imin. 



Wl = READPRN (link l_wgt) 
W2 = READPRN (link 2_wgt) 
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9x(5+l) weitht matrix between input 
layer and the first hidden layer 



W2 = (0.43 1.46 1.39 3.91 3.96 2.16 2.25 -6.36 -4.1) 2x9 weight matrix between 

hidden layer and output 
layer 
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